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ABSTRACT 

Context. Starburst galaxies show a direct correlation between radio and far-infrared emission. High target densities and a high rate 
of supernova explosions imply the possibility of accelerating hadronic cosmic rays and producing decay products from hadronic 
interactions, like high-energy neutrinos and photons. 

Aims. We propose an explanation for the far-infrared/radio correlation of galaxies in terms of the energy balance of the interstellar 
medium and determine the flux from high-energy photons and neutrinos from starburst galaxies. 

Methods. We present a catalog of the 127 brightest starburst galaxies with redshifts of z < 0.03. In order to investigate the correlation 
between radio- and far-infrared emission, we apply the leaky box approximation. Further, we derive photon- and neutrino spectra 
from proton-proton interactions in supernova remnants (SNRs). Here, we assume that a fraction of the SNR's energy is transferred to 
the acceleration of cosmic rays. We also investigate the possibility of detecting Gamma Ray Bursts from nearby starburst galaxies, 
using the catalog defined here. 

Results. We show that the radio emission is only weakly dependent on the magnetic field. It turns out that the intensity of the radio 
signal is directly proportional to the number of supernova explosions, which scales with the far-infrared luminosity. In addition, 
we find that high-energy photons from proton-proton interactions in SNRs in starbursts can make up several percent of the diffuse 
gamma-ray background. The neutrino flux from the same sources has a maximum energy of ~ 10' GeV. Neutrinos can, on the other 
hand, can be observed if a Gamma Ray Burst happens in a nearby starburst. About 0.03 GRBs per year are expected to occur in the 
entire catalog. The true number is expected to be even higher, since we only include the brightest sources. The number of events per 
burst in IceCube varies between about one event and more than 1000 events. This provides good prospects for IceCube to detect a 
significant event, since the background for a GRB search is close to zero. 

Key words. Galaxies: starburst - Infrared: galaxies - Radio continuum: galaxies - Catalogs - cosmic rays - Neutrinos 



1. Introduction 

• Radio emission from galaxies is usually dominated by syn- 
chrotron emission from a population of non-thermal, energetic 
electrons in a magnetic field which permeates most of the in- 
terstellar medium. This radio emission is often spatially struc- 
tured, such as in the starburst galaxy M82, showing individual 
compact sources, which can be interpreted as fairly young su- 
pernova remnants ( Kronberg et al.L 119851: iKronberg & SrameH 
119851: iBartel et all Il987l) . The origin of these energetic elec- 
trons, a part of the cosmic rays, is thus expected to b e the young 
supern ova remnants, see iBaade & Zwickvl (Il934l): IShklovskil 
(Il953h and for a extensive review Berezinskv et al.l ( fl990l) . Thus 
the radio emission is a key to interpret the physics of cosmic 
rays, and conversely, any attempt to understand cosmic rays 
should also try to understand the properties of the radio emis- 
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Galaxies also have abundant far-infrared (FIR) emission, which 
is due to dust. This dust is heated by stars, often mostly young 
stars. As was noted from the mid-eighties, this thermal dust 
emission correlates rather well with the non-thermal radio emis- 
sion. The correlation in its most simple form is just a propor- 
tionality between far-infrared and non-thermal emission. As ref- 
erence wavelengths, 60 //m and 100//m are used for the far- 
infrared. Frequencies between 1 .4 GHz and 5 GHz are used as 
typical for the radio regime. 

Many attempts have been made to understand the proportional- 
ity between radio and far-infrared. On the basis of rather simple 
modeling of galactic evolution, a strong correlation is actually 
expected, since both supernova remnants and the dominant heat- 
ing by ultraviolet l ight from massive stars derives from the sam e 
stellar population dBiermannl 119761 : iBiermann &FrickeL[T977h . 
Using such models t he far-infrared luminos i ty of NGC2146 
had been predicted by IKronberg & B iermann (1981) and veri- 
fied subsequently by IRAS observations (Moshiret al., 1990a]). 
The correlation as seen in the data was first clearly stated by 
Ide Jong et al.l (Il985h . and subsequently discussed at some length 
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by many authors jBicav et al.L [19891: IWunderlich et all 119871: 
IWunderlich & Kleinl [T988L 119911: ICondon et all [T99T1) . It still 
defies a clear explanation. 

An extensive attempt to int erpret the correlation was mad e 
by Volk and collab orators jVolkl 119891: IXu et all Il994bllal: 
lLisenfeldetaILll996h . The electrons are believed to lose all en- 
ergy in this model and therefore, the correlation is calorimet- 
ric. However, this would predict an actual steepening of the ra- 
dio emission between the reference frequencies, an effect which 
is not seen in the data. The solution of a spatial mixture of 
cutoffs would still allow for a locally loss-dominated scenario. 
However, first spatially resolved observations of M 33 indicate 
that the local star-forming re gions typically have flat spectra, 
too (iTabatabaei et al.L 12007 alibi) . It is therefore likely that star 
forming regions are generally injection dominated. 
Here, we propose a simple model, based upon a particular pic- 
ture of the energy balance in the interstellar medium. This model 
also uses some simple assumptions, as we will emphasize. The 
model is local, and so automatically allows for starbursts and 
gradients in disk galaxies, while upholding the correlation. The 
model leads to some specific predictions which can be checked 
with further data. A catalog of starburst galaxies is presented to 
perform first checks. 

The outline of this paper is as follows: We first define variables, 
used throughout the paper, in Section |2] In Section |3] a cata- 
log of 127 nearby, bright starburst galaxies is presented, includ- 
ing far-infrared, radio and X-ray data. The correlation between 
FIR and radio emission is outlined together with the difficulty 
of explaining it. In Section|4] we present a model explaining the 
FIR-radio correlation. The possible emission of cosmic rays and 
secondaries produced in proton-proton and proton-photon inter- 
actions is discussed in Section [5] In particular, we examine the 
possibility of detecting secondaries from supernova remnants, as 
well as cosmic rays and secondaries from Gamma Ray Bursts in 
starbursts. Finally, imphcations are discussed in Section |6] 

2. Definitions 

In the following sections, the electromagnetic spectra at different 
wavelengths are used in order to investigate the starburst nature 
of the catalog sources. Table [Ugives a summary of the different 
parameters used. Concerning spectral power-law fits between the 
wavelengths, we use the convention 
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with S as the flux per area and frequency interval, in units of 
Jy= 10"^^ W/m^/Hz. Here, Sq is the flux at a reference fre- 
quency Vq. 

3. Local starbursts: a sample 

In this section, we present a sample of local starburst galaxie^. 
The data of the individual sources are presented in appendix lAl 
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Fig. 1. 60fi luminosity - distance diagram. Crosses indicate all 
309 pre-selected starburst galaxies, squares show those remain- 
ing after the cuts S i.4ghz > 20 mJy and S eo/j > 4 Jy and z < 0.03. 
The dashed line shows the sensitivity for Seo^i > 4 Jy. 



Only local sources are considered, since our aim is to investigate 
the closest sources. The catalog presented here consists of a total 
of 127 starburst galaxies. This is a sub-sample from a larger sam- 
ple of 309 starburst galaxies, applying cuts at both FIR and radio 
wavelengths to ensure a complete, local sample. These cuts are 
discussed in the following paragraphs. Different tests were per- 
formed in order to verify that the considered galaxies are indeed 
starbursts, as presented in the following paragraphs. In order to 
remove contamination from Seyfert galaxies, we only use high 
ratios of FIR to radio flux density, i.e. 5'60/j/5' i.4ghz > 30. To test 
if the sample consists of starburst galaxies as opposed to regular 
galaxies, we check that the correlation between radio power and 
FIR luminosity is a direct proportionality. Apart from that, our 
main criterion for the catalog is that the sources are closer than 
z < 0.03, i.e. located in the supergalactic plane, and that they 
have both radio and IR detections. The latter gives us informa- 
tion about the ratio of the IR to radio signal, which we require to 
be larger than 30. This ensures a high IR component compared 
to the radio part, i.e. that the sources are indeed starbursts and 
not Seyfert galaxies. Further, we apply sensitivity cuts, we only 
include sources with a flux density > 4 Jy at 60 /jm and a radio 
flux density at 1 .4 GHz larger than 20 mJy. Figures[T]and[2]show 
the 60 yum resp. 1 .4 GHz luminosity of starburst galaxies versus 
their luminosity distance. The dotted lines represent the sensi- 
tivity for 4 Jy, resp. 20 mJy. Crosses represent all 309 sources 
we selected in the beginning, squares show those 127 sources 
remaining after the cuts at 5'i.4ghz > 20 mJy and Seo^ > 4 Jy, 
as well as z < 0.03. We apply those cuts in order to ensure a 
complete, local sample in both FIR and radio wavelengths. 
Since the sources are closer than z - 0.03, many of the star- 
bursts are located in the supergalactic plane. Their spatial dis- 
tribution should therefore be a flat cylinder with a further more 
spherical component, for those sources not in the supergalactic 
plane. We therefore expect that the number of sources with a 
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Table 1. Summary of parameters used in this paper. 
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Fig. 2. 1.4 GHz luminosity - distance diagram. Same notation Fig. 3. log A' - log 5 representation of the catalog. An S '^-fit 
as Fig. [1] The dashed line shows the sensitivity for S i 4ghz > matches the data nicely, with a turnover at 5o = 10.56 Jy. 
20 mJy. 



flux density larger than S , N{> S), should follow a behavior of 
- 5"''^. A pure 5 "'-behavior is expected for a flat cylin- 
der, while a spherical distribution results in an 5 ' ^ -behavior. 
Figure [3] shows the logarithmic number of sources above an FIR 
flux density Seo/i- We fit the data with the following function: 



N{>S)^No-(S +So) 



-P 



(2) 



Here, No, So and /3 are fit parameters. Using an error of "^V, the 
parameters are determined to 



No 
So 



3155 ± 1297.9 
(10.56 ± 3.78) Jy 
1.2 + 0.2. 



The behavior N{> S) ~ ^±0.2 [patches the expectation that 
the function should lie between 5 ' " and 5 In the following 
paragraphs, we will investigate further whether the classification 
of the 127 sources as starbursts is justified. 
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Radio Luminosity vs. FIR Luminosity (z>0.03, S60^m/S 
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Fig. 4. Radio power f i.4ghz at v = 1 .4 GHz versus FIR luminos- 
ity LpiR- A direct proportionality, Pi.4ghz ^ Lfm is found. 

3.1. FIR luminosity versus Radio power 

Looking at a well defined sample of galaxies, it turns out that the 
correlation between radio and far-infrared (FIR) emission is not 
linear, i.e., that the radio luminosity is prop ortional to the far - 
infrared luminosity to the power 1.30 + 0.03 jXu et aUri994bh . 
As Xu and collaborators note, the far-infrared emission has two 
heating sources, stars that do explode later as supernova rem- 
nants, and also stars, that will never explode as supernovae. This 
second population of stars needs to be corrected for, and their 
contribution to the dust heating needs to be eliminated. This 
then leads to a corrected far-infrared luminosity, which is di- 
rectly proportional to the radio luminosity (Xu et aL, 1994b). 
The proportionality holds along a disk in a galaxy, even for fairly 
short lived phases like a starburst, such as in M 82, and thus re- 
quires clearly local physics, with a short readjustment time scale. 
This poses a severe difficulty for any proposal to explain the ra- 
dio/FIR correlation. 

Th e FIR luminosity in the range of 60 /zm and lOOyum is given 
as (IXuet al.LfT994bl) : 

LpiR :=4;rfl'f -Ffir. (3) 
Here, c// is the luminosity distance of the individual sources and 



:= 1.26- 10" 



2.58- 



> lOO/j 
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is the FIR flux density at Earth as defined in lHelou et al.l (Il988l) . 
The normalization factor comes from the frequency integration 
and from the conversion of Jy to W/m-/Hz. In Fig. |4l the loga- 
rithm of the radio power at 1 .4 GHz, Pi .4ghz versus the logarithm 
of the FIR luminosity Lfir is shown for our catalog. The circles 
show the single sources and the solid line is a fit through the 
data. The fit yields a correlation of 
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This demonstrates that short stellar lifetimes dominate the cor- 
relation in our sample, and so this is strongly supporting our hy- 
pothesis, that the majority of our sample galaxies are starbursts. 

3.2. Infrared to radio flux density ratio 

Generally, regular galaxies are distinguished from active galax- 
ies by their ratio of the FIR flux density at 60;um, Sa)fi, and the 



radio flux density at 1 .4 GHz, S i.4ghz: 



■seo/i/ 1.4GHz 



>60/j 



' 1.4GHz 
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For Seyfert galaxies, this ratio is about sa)fx/iACiiz ~ 10, while 
it is significantly higher in the case of starburst galaxies, 
i60/j/ 1.4GHz ~ 300. The histogram of the ratio between the FIR 
flux density at 60yum and the radio flux density at 1.4 GHz is 
shown in Fig.|5] All 127 sources have a ratio of .s'60/i/i.4GHz > 30, 
which confirms that the sources are not likely to be Seyferts. 
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Fig. 5. Ratio of the flux density at 60 //m and at 1.4 GHz. AU 
sources in the sample have ratios larger than 30, which indi- 
cates a high star formation rate. The me dian is around 1 0. Thi s 
matches previous investigations, e.g. iBiermann et al. I (II985h . 
who find a mean value of 250 at higher radio frequencies, 
V = 5 GHz. 



3.3. Radio to Infrared and X-ray to Infrared spectral indices 

A further criterion of distinguishing regular galaxies and 
Seyferts is their spectral index from X-ray to IR (XIR) 
and from radio to IR (RIR). The diagram of the XIR 
(I keV to 60 jum) versus RIR (5 GHz to 60 jum) index of 
the sources in shown in Fig. |6] Derived from figure 3 in 
(iRodriguez-Pascual et al.L Il993l) . starburst galaxies have spec- 
tral indices scattering around (RIR, XIR).,,^^/:,!!^ ~ (0.6, -1.9), 
Seyfert-I galaxies show (RIR, XIR)^,.-/ ~ (0.48, -1.2), Seyfert- 
II galaxies have (RIR, XIR)^,,-// ~ (0.47, -1.6) and quasars 
are located at (RIR, XlR)guasar ~ (0.28, -1.1). The values for 
the RIR and XI R indices of starburst galaxies given by 
IChini et al] ( Il989[) a re slightly h i gher, w hich matches the sam- 
ple examined here: IChini et akl ( Il989l) give a RIR index of 
0.82 and a XIR index of -1.66. We find average values of 
(RIR, XIR) = (0.82, -1.77) which is compatible with the ex- 
pected result. 

Still, we do not have X-ray data for all the sources, so there 
may still be some contamination from both Seyferts and regu- 
lar galaxies in the sample. As we only used catalogs where the 
sources have previously been identified as starbursts, this con- 
tamination should be small. 
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Fig. 6. Radio-to-IR spectral index versus X-ray-to-IR spectral in- 
dex. The crosses represent those 48 sources in our catalog with 
radio, FIR and X-ray measurements. The blue triangle shows 
the average of the values. The open circle shows the average lo- 
cation of Seyfert-I galaxies, the open square represents average 
Seyfert-II galaxies a nd the star indicates QSOs. The last three 
values are taken from lRodriguez-Pascual et aTl ( |1993[). Note that 
indivi dual galaxies scatter around the given values dChini et al.L 
119891) . 



4. The interstellar medium and the FIR/radio 
correlation 

The interstellar medium connects the formation of stars, the 
explosion of supernova remnants, the regularization and en- 
hancement of the magnetic field, and the transport of cosmic 
rays. In order to understand the observation that the thermal 
hot dust emission from a galaxy is simply proportional to the 
non-thermal radio emission from relativistic cosmic ray elec- 
trons, we need to understand the interstellar medium, or at least 
get close enough to a comprehension, that we can understand 
this amazing correla tion (dejong et al., 1985; Wunderlich^et al.„ 
[19871 1 WunderUch & Klefnl i 1 98i i l 99 it fCondon et all 1199 iT) . 
Obviously, since the very massive stars power the far-infrared 
emission through a large fraction of the ultra-violet emission, 
which is absorbed by dust, and then in supernova explosions pro- 
duce the energetic cosmic ray electrons, which emit the observed 
radio emission, there should be a correlation. Using simple ini- 
tial models for stellar population evolution, this rather naive 
early picture successfully predicted the app roximate far-infrared 
emission already in 197 7 ('Biermann'. ll976tlBiermann & Frick5 , 
Il977t iKronberg & Biermann. ,198L) . However, the tightness of 
the correlation was neither anticip ated nor predi c ted. T he cor- 
relation was finally discovered by Ide Jong et al.l ( Il985h . Since 
the dust emission just measures the total output in ultra-violet 
by young stars, it seems obvious, that in the limit of much ab- 
sorption, the far-infrared emission would just be proportional to 
the star formation rate, given a general initial mass function. 
However, the radio emission is approximately proportional to 
the product of the cosmic ray electron density and the magnetic 
field energy density, and therefore it is not really obvious at all, 
that integrating along a vertical column through the disk of a 
galaxy these two emission components should be basically pro- 
portional. 

Modern descriptio ns of starbu r st galaxies su ch as M 82 or 
NGC 2 146 are in iDopita et al.l (I2005L l2006allbh : iGroveset al.1 
(l2008h . 



It had been observed early that the three main components of the 
interstellar medium, the gas, the cosmic rays, and the magnetic 
field have very similar energy densities, or pressures. Since all 
three derive from very different physical processes, to keep them 
at approximate equipartition implies that the three time scales of 
change are also all three the same. This is the basic premise of 
the following argument, and it will lead naturally to an under- 
standing of the far-infrared radio correlation. So in t his spe cific 
sense it is a calorimetric argument similar to IVolkl d 1 9891) . al- 
though we approach the problem is a somewhat different way. 
The following line of reason i ng is v isibly influenced by 
earlier 
) 



like 



1951 



1977); 



paper s 
; ICoxl jI972h 



Biermann (1950); Biermann & S chliitei 

lCox & Smith ( 1974) ; McKee & Ostrikei 

Beuermann et all i 19851) ; 



■Kronberg et al.l 
(Il985h; ISnowden et all 




Kronberg & Sramek , 

Hunter etal. (199^; iBeck et alJ (l2003h ; iHanasz et al.l 
2006), and of course others mentioned in due course. 
First we wish to establish the concepts which we use, for easy 
reference, and then apply them to the problem here. 



4. 1 . The three main components of the interstellar medium 

The three main components are the gas, the cosmic rays and the 
magnetic field. At least the gas and the magnetic fi eld is clearly 
spatially highly inhomogeneous dBeck et al Ll2003h : 
The gas has a number of components, molecular clouds, neu- 
tral Hydrogen clouds, diffuse neutral Hydrogen, diffuse ionized 
Hydrogen, HII regions, stellar wind bubbles, supernova rem- 
nants with X-ray emitting shells, a tunnel netw ork of connected 



older supernova remnants dCox & Smith , 



( Beuermann et all 1 1 9851: Snowden et al 



120071) . and a wind dBreitschwerdtl l2008i lEverett et al.L \20m 



1974 ) a thick hot disk 
1997t IKronberg et all 



The tunnel network probably connects to the hot thick disk. The 
wind is probably fed from the hottest regions of the tunnel net- 
work, in a fashion perhaps similar t o the Solar wind being fed 
from coronal holes (IStepanian et all [20081 e.g.). 
The magnetic field is permeating almost everything, and also has 
a thick disk. The field permeates the clouds, and is often con- 
fined visibly by the clouds (Appenzeller, 1974, e.g.). The mag- 
netic field is strongly perturbed by HII regions, and supernova 
explosions. The magnetic field is transported out of the disk by 
the wind. 

The cos mic rays, produce d by su pernov a explo s ions in shock- 
waves dBaade & Zwickvl 1 1934 iFermI [19491; [Prurvl [1983[; 
iBerezinskv et all 1 19901) cannot easily be repelled by anything, 
and so go through all clouds, and all the neutral and ionized gas. 
In bulk they cannot travel faster then the Alfven speed, since 
otherwise they would excite waves in the plasma, scattering the 
particles, effectively reducing their bulk velocity. 
The magnetic field is instrumental to allow cosmic ray accelera- 
tion in shocks, and perhaps throughout the medium. 
The cosmic rays in turn drive the dynamo mechanism to en- 
hance an exis t ing magnetic field and give it spatial coherence 
dParkeii [T969[ [1992); [Ferrierel [19961: [perriere & Schmii. [20001; 



iHanasz et al r [2004 [20061; lOtmianowska -Mazur et alllaOOSl). hi 
the classical Biermann-battery mechanism dBiermannl [19501 ; 
Biermann & Schluter, 1951), only a rotating star with surfaces 
of density and pressure in non-coincidence is required to pro- 
duce a seed field, which, however, is generally weak; the dynamo 
mechanism in stars can st rongly enhance magnetic field s, and 
through winds eject them dBisnovatvi-Kogan et all [r973L e.g.): 
This would constitute a very irregular, but potentially relatively 
strong source of magnetic fields in galaxies; in such a case the 
dynamo mechanism on a Galactic scale is required more to reg- 
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ularize the field rather than to strengthe n it. In the mechanism of 
iLucek & Belll(l2000h : lBeil (12004 12005h the shock waves can di- 
rectly enhance the magnetic field, using an existing population of 
cosmic rays. As noted already, the cosmic rays couple eff'ectively 
to the gas. Anisotropics in the phase space distribution of parti- 
cles, always present in Shock waves, a lso can produce new mag- 
netic fields on small scales ( IWeibej 119591 : iBvkov & Toptvgiru 
2005h. Galactic magnetic fields h ave also been reviewed b y 



Beck et al. (ll996h : lKulsrudl(ll999l) : lKulsrud & Zweibeil(l2008l) . 



Thermodynamically both the cosmic rays and the magnetic field 
can be thought of as a relativistic gas, with almost zero net mass 
density. 

Therefore the ensemble of cosmic rays and magnetic field con- 
stitute a light fluid pushing against the heavy fluid of the normal 
gas, and given enough energy density, these two components e s- 
cape via an instabihty ( |Parkeilll965tlKowal et all.l2003Ll2006l) . 

Therefore all three components are strongly coupled, and the 
data confirm an approximate energy equipartition between cos- 
mic rays and the magnetic field, and the sum of these two com- 
ponents equal in energy density to the gas. 

The one given parameter is the total energy input, integrated 
across all spatial inhomogeneities, since the energy supply is 
given by the stars, in the form of winds, explosions, and radi- 
ation. Since the inhomogeneities are extreme, especially in the 
density, it is important to use spatially integrated energy densi- 
ties for reference as much as possible. So we note that the energy 
density of the magnetic irregularities inte grated over all spa tial 
scales is also in approximate equipartition dBeck et anil996l) . 

The energy input can be estimated from the explosions of su- 
pernovae to about 1 supernova of 10^' erg, every 100 years, so 
at Lkin = 3 ■ 10^' erg/s; the uncertainty in this is about a factor 
of 3. Some fraction of this energy g oes into cosmic rays. This 
fraction could be large (lDrurvlll983l) . 

Othe r energy input can be estimated from the infrared emis- 
sion ( ICox & Mezgeii Il989h : About 1/3 to 1/4 of the total stel- 
lar radiation is absorbed by dust and reradiated in the infrared, 
beyond a wavelength of 25 ju, about 10'° L©. Of this, about 
2 ■ 10^ Lq = 8 • lO'*" erg/s is coming from young star forming re- 
gions. 

The energy density of magn etic fields can be estimated to 
be ab out 1.6 ■ 10 dyn/cm^ dBeck et al.L 119961: lEverett et all 
l2008h . flie energy density of cosmic rays is about the same, and 
their sum is about equal to the gas pressure, of 4 ■ 10"'^ dyn/cm^. 
Using a scale height of full width of 3 kpc, and a radius of 10 
kpc, we obtain a crude estimate of the energy content. This re- 
quires for magnetic fields and cosmic rays together an average 
supply of energy of about 4 • 10"*' erg/s, using the time scale ob- 
tained from cosmic rays (see below). Su ch numbers are uncer- 
tain by probably a factor of 2. Everett et a l. (2008) also estimate 
the required wind-power to about 4 ■ lO'*' erg/s, a Galactic wind 
driven by cosmic rays; they discuss other estimates. 

It is interesting to note that to within the uncertainties all these 
power estimates (supernovae, wind power, magnetic field and 
cosmic ray replenishment) agree better than their respective er- 
ror estimates. 

Therefore the time scales to replenish anyone of the components 
must also be approximately be the same. We do have the real 
number from radioactive isotopes of c osmic rays interacting, and 
the number is about 10 million years dBrunetti & Coding l2000l 
e.g.). This is then the time scale for all key processes. 



4.2. Supernova explosions 

For didactic simplic ity we basically adopt the approach of lSedovl 
dT958l) : ICoxl([T972l) . but use more modern cooling approxima- 
tions, and allow for much lower environmental densities, but 
otherwise rescale their equations. 

An explosion runs into the interstellar medium, and expands into 
the tenuous gas, which surrounds the clouds, and extends far 
above and below the central layer of cool clouds. We consider 
the expansion into the surrounding low density medium and ask, 
when the expansion runs into the cooling limit: 
The first question is what de nsity should be used: 



The galaxy has a wind ( Westmeier et aU 2005 



'Breitschwerdt', "200^ 
2008; Otmianowska-Mazur et al. 



lEverett et all 120081: iGressel et al 



^. , ^, , 20081) and it is getting fed 

from the tunnel network of ICox & SmithI d 19741) probably. This 
impUes that the Alfven speed must approach the escape speed, 
for a cosmic ray driven magnetic wind. Ot h er galaxies also 
show evidence for winds (C hyzv et al.Ll200()allbl : IChvzv & BecH 
,20041: IChvzv et al.L ,2006. 20^ : 



Va ^ — ^ 400km/s. 
4-np 



(7) 



With Z? a: 3 yuGauss, this implies a density of about 

« = 3 ■ lO-^cm-l 

The time to start cooling is 



= 5 ■ 10''yrs 



n-3.5 



2/11 



(rA_2i n-3.5) 



-5/11 



(8) 



where £51 is the energy of the explosion in units of 10^' erg, 
n-3 5 is the tenuous density in units of 3 ■ lO"'* cm"-', A_2i is 
the cooling coefficient in units of 10"^' erg cm"* s"', and r is a 
compaction parameter of order unity. This low density reflects 
the finding that the tenuous medium is of very low density, and 
due to substructure may on volume ave rage be of even lower 
density than s ugges ted by the X-ray data dSnowden et all 119971: 
lEverett et al.L l2008l) . of order 3 ■ 10"' cm"^; but we do use the 
temperature o f 10^ K, nea r the maximum, and also close to the 
stable region dFieldlll965l) . However, the cooling suggested by 
the X-ray spectrum is an integral over the entire evolution, and 
so only sensitive to the earliest part of the evolution. As soon 
as two or more supernova remn ants overlap, and start building 
a network dCox & SmithI 1 19741) . then the temperature evolution 
will be diff'erent, giving again higher temperatures, consistent 
with observations. 
The radius at that stage is 



i?'' = 8 ■ 10^ pc 



«-3.5 



3/11 



(rA_2i n-3.5) 



-2/11 



and the temperature then is initially 



T'^ = 1.7- 10^ K 



£51 
n-3.5 



2/11 



(rA-2i n-3.5) 



6/11 



(9) 



(10) 



This corresponds to an injection scale of turbulence. 
Interestingly, the time scale is of the same order of magni- 
tude to what we derive from cosmic ray transport, and the 
length scale is not fa r from the scale height of the hot disk 
dSnowden et all 119 97). demonstrating qualitative consistency. 
We note that lSnowd en et al. ( 1997) gave a much higher temper- 
ature, of about 4-1 0^ K, with a density of 3 ■ 10"' cm"^. Also, 
lEverett et al.l (l2008l) suggest a higher temperature. However, the 



Julia K. Becker et al.: Cosmic Rays VI 



7 



luminosity of that phase is a very small fraction of the entire 
dissipation in the ISM. 

As [Field (1965) shows the cooling is stable if the temperature 
dependence of the cooling function A is sufficiently strong and 
its double logarithmic derivative positive, or in the presence of 
heating larger than 2. This is the cooling phase we are inter- 
ested in. This is the case at temperatures below and near about 
~ 10^ K. 

Data suggest that the break-up of supernov a remnant shells ha s 
been observed in the starburst galaxy M82 dBartel et al.Lll987l) . 
However, in that case it is not clear whether we are observing 
the break-up of a wind-shell produced in a snow-plow effect by 
the stellar wind prior to the supernova explosion, or the break- 
up of the snow-plow of the normal supernova exploding into the 
interstellar medium. 



4.3. Magnetic inhomogeneities 

In the magnetic field data in our galaxy (iBeck et al.Ll2003l) there 
is already strong evidence for small scale substructure, since dif- 
ferent measures of the magnetic field yield very different num- 
bers: linear measures such as Faraday Rotation Measures indi- 
cate much lower strengths of the magnetic field than quadratic 
measures such as s ynchrotron emission. This is typical for small 
scale substructure jLee et al. Avillez & Breitschwerdd 

l200ll2007l) . where for a given total energy content high inten- 
sity sheets can hold all the energy for a small volume fraction. In 
such a picture linear measures give a much smaller number than 
quadratic measures, as is well known from mathematically iso- 
morphic arguments in thermal emission. Of course we should be 
comparing the proper integrals, also involving the spatial distri- 
bution of ther mal electron density a nd cosmic ray electron den- 
sity (see, e.g., iBowver et all (Il995h '). We ignore all this in our 
simple exercise. 

We can quantify this by integrating along a long thin cylinder of 
unit length. We refer to the magnetic field as Bq, when it is ho- 
mogeneous, and for the inhomogeneous case the magnetic field 
is Bi over most of the length, and enhanced by a factor 1/x in 
a region of length x: This then gives for the integrated energy 
density 



Bl--+B?-(l-x) = Bl. 



(11) 



Keeping the integrated energy content B^ constant, the linear 
measure of the magnetic field is given by 



Bl ■ - ■ X + Bl ■ (I - x) = Bl ■ (2 - jc) . 

X 

Combining the two expressions gives 



(2-x) 



(12) 



(13) 



for the ratio of linear measure versus quadratic measure. In the 
limit of small x this is just -y/I. The ob servations suggest that this 
ratio is of order 1/5 (iBecket al.ll2003h . and so x = 0.04 by order 
of magnitude. This implies that most of the magnetic energy is 
contained in shells of a volume a few percent, possibly as low as 
1 percent. Since the linear measure is proportional to the bending 
of ultra high energy cosmic rays, this implies that the bending 
is reduced by a factor between 5 and 10 over what we might 
reasonably expect otherwise. Obviously, in realistic situations 
much of this effect will be smoothed out, and so perhaps even 
more extreme situations may be required. 



Using the approach of ICoxl d 19721) with the env ironment of the 
tenuous hot phase of the interstellar medium (ISnowden et all 
1997; Everett et all l2008l) the cooling stage of an expanding 
shell of a supernova remnant might lead to such a configuration, 
of a very thin shell at large distances, with strong magnetic fields. 
In such a picture this stage would encompass most of the super- 
nova's energy dissipation, and so similar considerations may ap- 
ply to the interpreta tion of the X-ray data (iSnowden et al., 1997t 
lEverettet al.ll2008l) . 



4.4. Turbulence 

Turbulence is an ubiquitous phenomenon, and also i s a key 
ingred i ent in the interstellar medium (see reviews by iRicketj 
(1 19771) : [Goldstein et al.'('l995')). Key concepts to turbulence the- 
ory ha v e been introduced by Prandtl ( 1925); Karman & Howarth 



1938); 



1948); 



KoknogoroyI (Il941ailbl ig): lObukhovi ( Il94lh : iHeisenberg 
KraichnanT i 1965 ). and have been reviewed by Sagdeev 



19791) . One key argument which we wish to use, is the concept 



of the turbulent cascade. There the energy of the turbulence is in- 
jected into the gas at some large wavelength, and cascades down 
through wavenumber space, to the small wavelengths where the 
energy is dissipated. In many examples this leads in a three- 
dimensional isotropic model to the Kolmogorov cascade, which 
can be described in a local approxim ation by the fol l owing diffu- 
sion e quation in wavenumber space (iMcIvoii Il977t lAchterbeT^, 
Il979h : 



Here I{k) is the energy density of the turbulence per wavenumber 
k, and per volume element, and is the time scale of diffiision, 
which can be written as 



1 



Tk 



k{r,ffmk/p) 



1/2- 



(15) 



Here p is the matter density, and y^.//- is an effective adiabatic 
constant for the turbulent energy. The turbulence has a source- 
term, here limited to a single wavenumber fco. The turbulence dif- 
fusion equation basically says that the turbulence moves through 
wavenumber space with no additional source or s ink, as a con- 
stant energy current in wavenumber phase space dKolmogorovl 
J[941ai.baCj). The solutions to this diffusion equation can be writ- 
ten as 



m 
m 



k for k < ko and 
r^^^ for k>k„. 



(16) 



This latter behavior is commonly referred to as the Kolmogorov 
cascade, and is found ubiquitously in nature. 

4.5. The cooling of the interstellar medium 

The interstellar medium has a number of phases, which appear 
to be in approximate pressure equilibrium. Concentrating on the 
phase of the highest temperature, the highest speed of signalling 
(be it sound waves, or Alfven waves, or other wave modes), we 
note that its temperature is in the range where the typical gaseous 
emission is detected in the X-ray regime. The cooli ng curve of 
such a gas has been extensively discussed by many dCoxL [19721 : 
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[Sutherland & Dopital [T99llDopita & Sutherland! I2003L e.g.). It 
shows the following features, starting at low temperature, and 
considering the cooling coefficient A(T) with « the interstellar 
density in particles cm"-' 



n" A{T) erg /cm /sec . 



(17) 



This cooling curve A(r) rises from near 10"* K to a lo- 
cal peak near » 3 ■ 10^ K. The level of cooling along this 
local peak is given by A a; 10"^' ergs/cm^-'/s, dropping to 
A X 10"^^ ergs/cm^-'/s near x 10^ K, and towards a mini- 
mum near A ^ 3 • lO"^-' erg/cm^-'/s in the range x 3 ■ 10*" K 
to a: 10^ K. The peak is due to many edges and lines in 
the soft X-ray range. At higher temperatures the continuum 
emission begins to dominate and that emission is then given 
by A ^ 1.4- 10"^' r'^^ erg/cm^-'/s. The sharp cutoff to low 
temperatures near 10"* K is due to beginning recombination 
and thus a smaller density of free electrons to interact with. 
A contemporary discussion including the effects of strong 
departures from ionization equ il ibrium has been given by 
ISchmutzler & Tscharnuted (ll993h : lBreitschwerdt & Schmutzlej 
(11994 . 

4.6. The leaky box approximation 

Consider a column perpendicular to a galactic disk of height H, 
and the number of cosmic ray particles in it as a function of par- 
ticle energy and time A^(£', t); then we have the balance equation 



d N{E, t) 

at t(E) 



Q(E). 



(18) 



where t(£') is the escape time scale. The sign of the second term 
is positive, since the process described is a loss. In a stationary 
state we then obtain readily 



N{E) = Q{E)t{E). 
The characteristic time of loss can be written as 

riE) ~ ^ . 

K 



(19) 



(20) 



Here, k is the diffusion coefficient of cosmic ray particles, which 
can be written approximately in the quasi-linear approximation 
as 

'- ''^'^ (21) 



K ~ — r„c 



mk ' 



where we write for the turbulent energy density 

I(k)k ~ Uo{k/k,)-^" ~ U„B-^"rf''E^'\ (22) 

using here the assumption of a Kolmogorov spectrum. Here, rg ~ 
l/k ~ E/B is the Larmor radius of a particle of energy E under 
consideration, which gyrates in a magnetic field of strength B. 
The basic radius ~ 1 /kg corresponds to the injection scale of 
the turbulence. This then leads to a dependence of the diffusion 
coefficient on the various parameters 



£1/3^5/3^2/3^-1^-1 



(23) 



Therefore, we adopt the point of view that the leakage time scale 
t(E) is proportional to (relativistic) energy E^^^^, and so that the 
equilibrium density of energetic particles A^(£') is proportional 
to E^^^^ as well. 

There is a difficulty with this argument, which can be solved: 
The secondary to primary ratio of cosmic ray nuclei such as 



the ratio Boron to Carbon already give information as to the 
energy dependence of the leakage time, and such an analy- 
sis gi ves an energy dependence as £'-" ^±0.1 ( ]£ngelmann et al 



. dEnge 

19901 e.g). We have argued elsew here already (iBiermamu 11995 



Wiebel-Sooth et al.L [19951 1 19981) . that this reflects the energy de- 
pendence of the amount of target material seen for spallation. 
Most of the target interaction with heavy nuclei among the cos- 
mic rays happens near the most massive stars, the Wolf-Rayet 
stars, with the turbulence excited by the cosmic rays themselves, 
giving an energy dependen ce of the spallation secondaries o f 
£-5/9 (|Biermann[ |T998; Bier mann et al.Ll200ll:lBiermannil2006l) . 
consistent with the data, which give 

£■-0.54 (|ptuskinlll999l) . On 
the other hand, most of the gamma emission from 7r-zero decay 
arises from the interaction among the more numerous supergiant 
stars, the red supergiants, for which we suggest that the turbu- 
lence arises from instabilities. 

4. 7. The radio emission 

In order to calculate the radio emission from the total number 
of relativistic electrons in a column we first have to proceed to 
work out the cosmic ray loss time for electrons, second the cool- 
ing time for the tenuous hot medium, put them equal, and third 
calculate the radio emission per supernova event from a column 
in the disk. 

The maximum energy in the elemental distribution and spectrum 
of the cosmic rays is in protons near their rest mass energy, since 
their spectrum is steeper than 2 in energy. This means that the 
energy which has to be used in the expression for the leakage 
time is fixed. So, the leakage time is given by 



H^Inkn 



TCR 



BVirl'^E'/^' 



(24) 



The maximum energy content of the population in the protron 
spectrum is here near E - I GeV, so that 



tcr(E = 1 GeV) 



H^Ipko 



(25) 



in this case. 

The cooUng time is given by 



nkgT 
Tl 1 



1 



{B^I'&nllkBT^f MT*) 
1 



(26) 



in the approximation that A(r*)/r; is a constant, and assuming 
equipartition again. This is reasonable in the dissipation stage of 
supernova remnants when A approaches a double-logarithmic 
derivative of 2, below a temperature of 10^ K. Also, in that tem- 
perature range the cooling time is a minimum, for a given energy 
density. We noted above that the dissipation stage reaches those 
temperatures. 

Putting the two time scales equal then yields the relation 



1 



(27) 



Now, we wish to consider electrons, which emit radio emission 
at a certain frequency v, which gives the condition that 



1/2^-1/2 



(28) 
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and so, making use of the frequency and magnetic field depen- 
dence of Equ. ( [26] l and Equ. ( |28] ), the diffusion time for electrons 
which emit at v is given by 

Tc«, ~ H^E-'^'B-'I\-^I'I^K ~ B-''l'y-'^' (29) 
The em ission is from a spec t rum o f electrons of £'-2.42±o.04 ^j. 
jection (iBiermann &StromL[T993l) . and so the total radio emis- 
sion S y per supernova event in a column in a disk can be written 
as 

r. nl.71±0.02 p -0.71±0.02 D- 11/6 -1/6 /om 

^ ^-0.12±0.02y-0.88±0.02 

for the case that all supernova remnants produce the same num- 
ber of relativistic electrons. Here, we use an integration along a 
vertical column to obtain the total number of energetic electrons 
per supernova explosion. This is reasonable since we are using 
the adiabatic phase of supernova remnant evolution, where the 
energy is conserved. This argument leads to such a weak depen- 
dence on the strength of the magnetic field, that the resulting 
offset is smaller than the errors in the data. 
Therefore the radio emission is only very weakly dependent on 
the magnetic field, independent of all other parameters, and is 
directly proportional to the number of supernovae per time in- 
terval, and so also to the luminosity of massive stars. 
The cosmic ray electron data show that above a few 10 of GeV 
the spectrum is a lready in the loss limit, so steeper by unity 
(lKardashevl[T962h . Below that energy the direct data are com- 
promised by the Solar wind modulation, but from the radio emis- 
sion of normal galaxies we do know that the electron spectrum is 
closely in agreement with what we obtain for protons, at some- 
what higher energy, allowing us to conclude that the electron 
spectrum corresponds to the diffusion limit, when the leakage 
out of the Galaxy is faster than the synchrotron loss. Using the 
expressions for the two time scales we can check: 
Radio emission at 5 GHz corresponds to an electron energy of 8 
GeV, and a synchrotron loss time of 4 ■ 10^ yr. 
As we note elsewhere, galaxies and starburst galaxies are some- 
times so young as to be in the injection limit, or so old as to be 
in the loss limit, rather than the leakage or diffusion limit that 
we consider here. For any reasonably young age of a starburst 
there is always some low radio frequency v\, below which we 
are still in the injection stage, below which the age of the star- 
burst is longer than the diffusion time scale. So below that the 
radio flux density is lower than the diffusion limit derived above, 
an d the radio flux dens i ty has the spectrum of injection, follow- 
ing lBiermann & StromI (Il993h 5 y ~ v , and so relative to the 
diffusion limit the radio emission at frequency v < vi is weaker 
by the factor (v/vi)'^^; this is typic ally not far b elow unity. At 
the the other extreme, considered bv lVolkl ( Il989l) there is also a 
radio frequency vo, beyond which the synchrotron and inverse 
Compton losses cut in, and become faster then the diffusion. 
Thi s is directly visible in the observed cosmic ray electron spec- 
tra jWiebel- Sooth & Biermannl Il999h . So, at radio frequencies 
V > V2 the radio emission is again weaker, this time by (v/v2)"''^. 
So again the radio emission is slightly down. Both variants show, 
that the radio emission is not far from the diffusive equilibrium. 
This concludes the demonstration of the argument. There are 
many checks on these ideas, which one can make, such as the 
pressure of the interstellar medium, the X-ray luminosity, the 
possibility to account for extreme galaxies s uch as M82 , the ra - 
dial gradient of the far-infrared/radio ratio ( iBicav et al.L Il989h . 
the thickness of the hot gaseous disk and the associated diffu- 
sion coefficient of cosmic rays, and many others. We will discuss 
these points elsewhere. 



4.8. Where is the limit to the diffusion iimit? 

The line of reasoning seems to be a deus ex machina in the sense 
that we seem to be always in the diffusion limit, independent of 
whether we have starburst galaxies with relatively strong mag- 
netic fields, or normal galaxies like our own. 
The equations above imply that 

= -^gi < - ^syn. (31) 

We note first that both sides depend on the magnetic field 
strength squared, so the argument is independent on the mag- 
netic field; Inserting numbers we find 

^ < i-i^, (32) 

A-21 Je 

where is the Lorentz factor of the mainly emitting electrons. 
This time, however, we need to ask what density contrast the 
main radio emitting substructures have, since the left hand side 
of the equation uses grand total volume averages, and the right 
hand samples just those regions, where B~ is especially high. 
Already simple arguments, as shown above, demonstrate that the 
average of is quite a bit higher than the average of B squared. 
We can assume here, that such a factor might be of order 30, or 
even higher, bringing the limiting energy of the electrons down 
to Je < 10'', or possibly even much less. We have already argued 
earlier that for the cooling to be maximal, the temperature has to 
be close to or below 10^ K. 

This allows us to understand perhaps, why galaxies are usually 
in the diffusion limit, sometimes of course, in a starburst, in the 
injection limit: In that case the overall spectrum corresponds to 
injection. There was just not enough time to achieve diffusive 
approximate equilibrium. 

4.9. Impiications 

First of all, this theory does just what one would naively expect, 
relate the massive stars which heat the dust through their ultravi- 
olet light directly with the subsequent supernovae. All the theory 
does, is work out the non-linearity inherent in synchrotron emis- 
sion, and shows them to introduce negligible dependencies on 
various parameters. 

Starburst galaxies as well as quiescent galaxies equally obey 
the correlation between the radio and the far-infrared emission. 
Therefore the theory implies by necessity that the magnetic field 
in a starburst region rises with the overall energy density. 
Galaxies which are subject to substantial compression by an en- 
counter with another galaxy, clearly have a magnetic field which 
is higher than corresponds to the energy density derived from 
star formation, and so may be expected to have more radio emis- 
sion than indicated by the general radio-far-in frared correlation . 
This is borne out by at least one example dHummel & Beckl 
[T995h . 

Furthermore, it is clear that for very short time scales, near order 
10^ years, the correlation cannot hold, since we then run into the 
lifetimes of the massive stars, which drive the energy balance in 
the interstellar medium. 

There is likely also a lowest level of star formation activity, 
where the assumption that the hot medium is fully connected, 
fails. There one may expect also substantial departures from the 
correlation. 
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5. Cosmic rays and their secondaries 

As stated above, the cosmic ray intensity from starburst galax- 
ies scales with the radio and infrared emission of the sources. 
In this section, we discuss the emission scenarios for charged 
cosmic rays and hadronic interactions leading to high-energy 
photon and neutrino emission. There are two source classes 
within starbursts that can accelerate cosmic rays to high en- 
ergies, namely shock fronts of supernova remnants and long 
Gamma Ray Bursts (GRBs), the latter being connected to su- 
pernova Ic explosions. In the first case, maximum energies are 
limited to less than 10'^ eV and thus, the cosmic rays from star- 
bursts cannot be observed directly due to the high cosmic ray 
background in our own Galaxy. Gamma Ray Bursts, on the other 
hand, were proposed as the ori gin of cosmic r ays above the an - 
kle, i.e. £cR > 3 ■ 10"^ eV, see lVietril (Il995l) : IWaxmanI (Il995l) . 
Since a high star formation rate as it is present in starburst galax- 
ies, leads to a high rate supernova explosions, an enhanced rate 
of long GRBs is expected. Thus, for closeby sources, the distri- 
bution of starburst galaxies can be used to test t he hypothesis of 
cosmic rays from starbursts, as also discussed in lBiermann et al.l 
(120081) . 

The dominant source of secondary cosmic rays like high-energy 
photons and neutrinos is proton-proton interactions in dense hy- 
drogen regions and proton-photon interactions in Gamma Ray 
Bursts. Proton-proton interactions produce pions via 



PP 



(33) 



Photohadronic interactions, on the other hand, yield pions via 
the Delta resonance. 



(34) 



High-energy photons and neutrinos are subsequently emitted in 
n"^- resp. ;7r°-decays: 



n yu -> e 



n 



yy ■ 



(35) 



5.1. Supernova remnants 



Cosmic Rays are believed to be produced in young super- 
nova remnants (SNRs), reaching maximum energies of around 
10'^ eV or above, depending on their local environment and on 
the cosmic ray composition. The production of secondaries from 
hadronic interactions depends on the proton-proton optical depth 
in the SNR environment. In this section, we present a model of 
which sources are optically thin and which, in contrast, are good 
candidates for the production of high-energy photons and neu- 
trinos. 



5.1.1. Optical depth 

We use the observation of synchrotron radiation from shock- 
accelerated electrons, assuming that electrons and hadrons are 
accelerated in the same shock environment. Figure [T] shows 
the distribution of spectral indices at radio wavelengths be- 
tween 1.4 GHz and 5 GHz for 105 sources in the sample with 
given spectral indices at the required wavelengths. The spectrum 
at these energies is produced by electron synchrotron losses. 
Depending on the shape of the primary electron spectrum and 
scattering effects, the spectral index of the electron population 



Spectral indices 5 GHz < V < 1 .4 GHz (z<^0.03. S60(Im/S >^30, S6()^lm>^4 Jy. S >^20mjy) 
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Fig. 7. Histogram of the radio spectral indices of 105 sources 
between 1 .4 GHz and 5 GHz. The areas between the dashed lines 
indicate sources in the loss, leakage and injection limit (from the 
left). As a conservative estimate, we include those three sources 
with extremely steep spectra as loss Umit sources. Those sources 
with extremely flat spectra, to the right of the injection limit area, 
are dominated by absorption or free-free radiation. 

can reach diff^erent values. Shock acceleration of charged parti- 
cles usually results in primary electron spectra of 



oc Ee 



(36) 



with aj'™ ^ 2.0 - 2.4. If the electrons escape before interacting 
with the ambient medium, the primary spectrum stays unmod- 
ified, a^"™ = af^, referred to as the injection limit. If the elec- 
trons are partly scattered down to lower energies, the spectrum 
of secondaries steepens to af^ « 2.5 - 2.8 for a primary spec- 
trum with ap'™ R! 2.0 - 2.4. This is called the leakage limit. In 
the case of calorimetric sources, basically the entire energy is 
lost in the source and the spectrum of secondary electrons is as 
steep as af^ x 3.2-3.4. This scenario is referred to as the loss 
limit. 

As discussed in iRvbicki & LightmanI (Il979l) . the electron spec- 
tral index af^ correlates with the index of synchrotron radiation 

as 



a - -- 



Thus, the observed synchrotron spectral indices are 



-0.75 in the injection limit 
— 1.0 in the leakage limit 
-1.2 in the loss limit . 




(37) 



(38) 



The synchrotron spectral indices in this sample scatter between 

-1.7 < Q- < 0.7 with a peak at ff 0.7. We exclude sources 

that may include contributions other than synchrotron radia- 
tion, i.e. those sources that have spectral indices a > -0.5. 
Here, absorption is likely to have modified the spectrum. 
Alternativ ely, spectral indices around -0.1 point to free-free 
radiation (iMezger & Hendersoni Il967h . Of the remaining 85 
sources, 37 sources (44%) starbursts are in the injection limit, 
36 sources (42%) are in the leakage limit and 12 sources (14%) 
are in the loss limit. 

The observation of a synchrotron spectrum following the in- 
jected electron spectrum has severe implications for hadrons in 
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the same source. As proton interactions require much higher par- 
ticle or electromagnetic field densities than electron interactions, 
it can be expected that hadrons are not affected if electrons can 
escape freely. Thus, protons escape from the source before they 
interact. Even in the leakage limit, it is unlikely to have pro- 
ton interactions with matter or photons, only sources in the loss 
limit provide conditions with reasonable densities for proton in- 
teractions. In the case of high particle densities, protons interact 
with each other and produce pions as described in Equations (l33T l 
and ( [35] l. Hence, only those 14% of all starburst sources in the 
loss limit are likely to produce both high-energy neutrinos and 
photons. We will take this into account in all following calcula- 
tions by applying a factor of e/oj, - 0.14, so that only loss hmit 
sources are included. 



5.1 .2. Contributions of starbursts to tine FIR background 

In a previous estimate of high-energy neutrino radiation from 
starburst galaxies by dLoeb & Waxmanl |2006|) . it is assumed 
that the entire background of far infra-red radiation come s from 
starburst galaxies. However, it is pointed out by iSteckeil ( l2007h 
that only a fraction of 23% the total diffuse FIR background 
actually originates from starbursts. This would reduce the flux 
of hadronic sec ondaries by a fa c tor ^f jR = 0.23. Still, as it is 
pointed out by [Thompson et al.l ( |2006|) . a fraction near unity, 
^FiR ~ 1, is consistent with most star formation rate models to- 
day and should be considered as an upper limit estimate. In this 
paper, w e will use the FIR b ackground from the EBL model 
given by dKneiske et al.Ll2002l) . where a total of 80% comes from 
starbursts. 



5.1.3. Production of liadronic secondaries 

The non-thermal radio emission from starburst galaxies indi- 
cates the shock acceleration of electrons. Hadrons are acceler- 
ated in the same way. In the case of proton interactions, high- 
energy neutrinos and photons can be produced as indicated in 
Equations ( |33T l and ( [35] l. For heavier nuclei, the efficiency is 
slig htly less than for pr o tons due to phot o -disintegration, see 
jHooper et all (l2005l) : lAve et al.l (l2005l) : lAnchordoqui et al.l 
In this paper, we calculate secondary spectra from pro- 
tons, and consider heavier cosmic rays elsewhere. High-energy 
photon emission can also be due to bremsstrahlung or in- 
verse Compton scattering, but above 1 GeV, pion dec a y pho - 
tons should dominate as discussed in iPaglione et al.l (Il996l) : 
iDomingo-Santamarfa & Torres! ( l2005b . Here, we re-calculate the 
possible neutrino and photon fluxes for those sources optically 
thick to proton-proton interactions. Normalization, spectral be- 
havior and propagation are treated as follows: 

- Spectral behavior 

We assume that the protons at the source follow a power- 
law spectrum with an index Up and with exponential cutoff' 

at £'m.TX , 



dNp 
dE„ 



exp 



(39) 



The normalization of the spectrum Ap is determined by as- 
suming that a fraction of the total SNR energy, 77, goes 
into cosmic rays. This is described in more detail in the 
next paragraph. The observed spectrum of cosmic rays be- 
low the knee, i.e. below 10'^ eV is typically assumed to 
come from supernova explosions in the Galaxy, with a cut- 
off at 10'^ eV. The spectral part between the knee and the 



ankle of the cosmic ray spectrum is still a matter of de- 
bate, but may arise from the heavy particle component from 
SNRs, see (Stanev et al., 1993). Here, we assume that SNRs 
in starburst galaxies produce similar spectra. Protons are 
likely to have an energy cutoff even below 10'^ eV. As the 
exact cutoff energy is not known, we use /Smax = 10'^ eV 
as an upper limit. The observed spectral index is Up - 
2.7. Stochastic particle acceleration usually produces spec- 
tra of around 2.0 < Up < 2.4. Hence, it is not sure yet where 
the steepening of the spectrum occurs, whether it is an 
internal steepening or a propagation effect. We therefore 
test proton spectra with indices of ap - 2.7, 2.4, 2.2, 2.0. 
A spectrum of E'^-'' would be present when diffusion ap- 
phes, spectral indices between 2.0 and 2.4 are predicted 
by stochastic acceleration without significant diffusio n. This 
approach differs from the calculations by Paglion e et al.l 
(Il996l) : iDomingo-Santamarfa & TorresI (l2005h . who use the 
diffusion-loss equation to determine the spectral index. 

In all following calculations, we determine the spectra of 
hadronic secondaries from proton-proton interactions, us- 
ing the delta-functional approximation for proton energies 
Ep < 100 GeV, see e.g. (Mannheim & SchUckeiser, 19'94- 
At higher en ergies, the more exa ct analytic approximation as 
presented in lKelner et al.l (l2006l) is used, where Monte-Carlo 
simulation results are approximated by analytical equations. 
For all energies, the logarithmic increase of the proton- 
proton cross section with energy is taken into account as 
described in iKelner et al.l (l2006l) . We have a high hydro- 
gen density, nn = 100 cm as we expect the dominant pro- 
ton acceleration and interaction to occur in heavy supernova 
remnants, having red supergiants or Wolf-Rayet stars as pro- 
genitors. 
Normalization 

We normalize the energy spectrum of cosmic rays by deter- 
mining their energy density pcR- The latter is defined by the 
energy integration over the differential cosmic ray flux, mul- 
tiplied with the energy, multiplied by a factor Anic to get 
from a flux to a density. 



An 



PcR : = 



dE„ 



dNp 
dE„ 



(40) 



In order to determine the normalization of the total proton 
flux from starbursts, we assume a power-law spectrum with 
an exponential cutoff as described in Equ. ( |39l ). 

The energy density, in turn, is directly proportional to the 
supernova rate in a galaxy, «sn, assuming that a fraction rj 
of the total energy of a supernova, E^nr ~ 10^' erg is trans- 
ferred to cosmic rays. 



PCR 



^loss ■ V ■ ^SNR • »SN 

c ■ driz) ■ {I + Z)^ 



(41) 



Here, we use eioss = 0.14 as discussed in Sections [5. 1.1 1 and 
15.1.21 We further assume that 5% of the total SNR energy 
is transferred to cosmic rays, i.e. = 0.05. To determine 
the supernova rate in a galaxy, the total FIR luminosity at a 
given redshift z, Lp]'j^(z), is estimated from the FIR emissivity 
Sfir{z), the latter comin g from the extragalac tic background 
fight (EBL) model from lkneiske et al.1 (|2002|) . 



dV, 



Si"'(z)—{z)dy. 
ctz 



Here, dV/dz is the comoving volume element. 



(42) 
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quantity 


variable 


value 


total energy release (SNR) 


EsNK 


10"' erg 


energy fraction transferred from SNR to CRs 


1 


0.05 


hydrogen density 


riH 


100 cm-^ 


fraction of starbursts in the loss limit 




0.14 



Table 2. Parameters used to determine the diffuse flux of high- 
energy photons and neutrinos from SNRs in starburst galaxies. 



The s upernova rate was determined by iMannucci et alj 
(l2003b to correlate with the FIR luminosity of the galaxy, 

«SN = (2.4 ± 0.1) ■ 10-'- ■ yr-i . (43) 

The FIR luminosity is expressed in terms of the solar 
luminosity, 

Lq - 3.839 ■ 10-^-^ erg / s. This relation was predicted 



in dBiermann & Frickel [T977HKronberg & Biermannll 19811) 
within a factor of 3, and other experimental results from 
Ivan Buren & Greenhouse! (1 19941) yield the same correlation 
within uncertainties. Finally, the cosmic ray energy density 
is given as 

PCR = 1.6- 10"'^ erg/cm^ ■ 

T] EsNR LfjR 



0.05 1051 erg 10^2 

diiz) 
6.5 Gpc 



(44) 



This number is compatible with the observed cosmic ray 
spectrum above 1 GeV in the Galaxy, - 1 eV/cm-^, as- 
suming a supernova rate of iisn = 0.03 /yr. 
- Propagation 

We apply that the particle energy at Earth is a factor of 1 + z 
higher than at the source. Neutrinos travel in straight lines 
without interaction. In the case of photon propagation, we 
include absorption e ffects by the EBL, using the model of 
iKneiske et alj ( |2002|) . 

From the considerations above, the spectra of hadronic secon- 
daries at Earth from a given redshift z are determined and a 
simple redshift integration is performed in order to get the to- 
tal neutrino flux at Earth. Here, we integrate from z = up to 
a redshift where first starbursts are expected to be formed. We 
use Zmax - 5, sincc the main contribution to the spectra comes 
from redshifts up to z ~ 1 - 2. The contribution above z = 5 is 
negligible. 



5.1.4. High-energy photons 

Assuming proton-proton interactions in starbursts as described 
above, we calculate the gamma-ray emission produced by star- 
burst galaxies for energies > 1 GeV. A detailed calculation 
and dis cussion for the starb u rst galaxy NGC253 has been 
done in iPagUone et al.l ( Il996l) : iDomingo-Santamarfa & TorresI 
(I2005h . Gamma-ray emission in starburst galaxies is due to three 
different processes, Bremsstrahlung, inverse Compton scatter- 
ing and pion decay. In the following only gamma-ray emis- 
sion from pion decay is considered be cause it is the domi- 
nating process at energies above 1 GeV ( iPaglione et aI.L 119961 : 
[Pomingo-Santamarfa & Torres, 2005). 

In Fig. [8] the model of gamma-ray emission from starbursts 
as presented above is shown for three different proton spectra. 



i.e. £-2 jjjg thick solid, £-2 2 ^jjg tjjiQjj^ dashed, £-2 4 
the thick dot-dashed and £-2 7 ^jje thick dotted line. We use 
that only 14% of all starbursts are in the loss limit {eioss = 0.14) 
and we apply the F IR background from the EBL model by 



apply the t n 

iKneiske et al.l (|2002|) . Absorption due to the EBL is consid- 
ered in all calculations. We compare our results to the obser- 
vation of the diffuse, extragalactic background as measured by 
EGRET ( Sreekumarj;t al., 199_8). Closed squares represent the 
first analysis presented in ISreekumar et al.l ( Il998[) . while open 
circles show an updated analysis, using an improved model for 
the distraction o f the galactic comp onent from the extragalac- 
tic contribution dStrong et al.L 120051) . The high-energy radiation 
from starburst galaxies makes up about 10% percent of the total 
background. In the same figure, we show other possible con- 
tributions to the total dif fuse backgro und, like the on e from 
resolved EGRET blazars (IKneiske &^ Mannh ei'ml l2008h or fli e 
contribut i on fro m regular galaxies ( Pavlidou & Fields! |2002|) . 
iKarlssonI (|2008|) prese nt a model of high -energy photon emission 
from starbursts. Just as lKarlssonI ( |2008|) . we include the logarith- 
mic rise of the proton-proton cross section with energy and get 
comparable results. In additi on, our results yield a sU ghtly lower 
flux than the prediction by iThompson et al.l (l2007h . The main 
reason is that we inclu de that only 14% of all s tarbursts are pro- 
ton calorimeters where pDiompson et al.l (l2007l) assume that this 
fraction increases from 10% for local starbursts to a saturated 
value of 80% at z = 1 . We use this more conservative approach. 
The considerations above show that several factors are quite un- 
certain, i.e. the total energy released by a supernova remnant, the 
hydrogen density, as well as the number of starburst galaxies that 
contribute to the FIR background and that are proton calorime- 
ters. Thus, the possibility of starbursts contributing significantly 
to the dilfuse photon background at high energies should be con- 
sidered. The Fermi satellitfl launched on June 11, 2008, will 
help identifying starburst galaxies at energies above 100 MeV, 
to determine the exact fraction of the starburst dilfuse photon 
flux in the total gamma background. 



5.1.5. High-energy neutrinos 

The high-energy neutrino flux calculated accord- 
ing to Section 15.1.31 is presented in Fig. |9] for an 
£-2", £-2-2^ £-2.4^ £-2.7 (jjjjgj^ gQjjj^ dashed, dot-dashed 

and dotted lines) spectral be havior of the initial pro ton spec- 
trum. Although, in contrast to lLoeb & WaxmanI (l2006h . we take 
into account that only 14% of all sources are effective proton 
calor imeters, the flux s trengt h of our prediction is coinpatible 
with iLoeb & WaxmanI ( |2006|) : While iLoeb & Waxmaiil ( |2006[) 
use the radio luminosity of local starbursts as an estimate, 
we take the supernova rate as a measure. The reason for the 
cutoff at lower energies is that we assume that protons are not 
accelerated beyond 10'^ eV, which is based on the observation 
of the knee in cosmic rays at about 10'^ eV. The interpretation 
that SNRs in the Galaxy are responsible for the flux below 
10'^ eV justifies the assumption that also in starbursts, SNR 
accelerate cosmic rays up to similar energies. It should be kept 
in mind, however, that the actual cutoff for protons is likely to 
lie at energies below 10'^ eV and that only heavier cosmic rays, 
which we do not consider in this calculation, can reach energies 
as high as 10''' eV. With an energy cutoff at or below 10'^ eV, 
the neutrino flux from starbursts is out of reach for the detection 
of a diffuse flux by high-energy neutrino detectors like Km3NeT 
and IceCube as indicated in Fig.|9l 



called GLAST before its launch 
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Fig. 8. Expected diffuse high-energy photon flux from SNRs in starburst galaxies. The black lines represent the contribution for an 
E^^ (thick, solid), E^^-^ (thick, dashed), E^^-* (thick, dot-dashed) and E'^-'' (thick, dotted) proton input spectrum with an exponential 
cutoff at 10''' eV. Closed data points are from the EGRET experiment (Sreekumaret aL, 1998). O pen d ata points represent an 
update of the EGRET data, using an updated mo del for galactic gamma - ray em ission (IStrong et al.L |2005|) . The thin, dashed line 
represents the contribution from E GRET blazars (Kneiske & Mannheimll2008l) . The thin, dot-dot dashed line shows the possible 
contribution from regular galaxies (iPavlidou & FieldsL I2002h r while the thin, solid lines display contribution from starbursts as 
calculated by [Thompson et al.l (|2007| ). E^^ and E^^-^ initial proton spectra. 
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Fig. 9. Expected diffuse high-energy neutrino flux from SNRs in starburst galaxies. The black lines represent the contribution for an 
E^^ (thick, solid), E^^-^ (thick, dashed), E^^-'^ (thick, dot-dashed) and E'^-'' (thick, dotted) proton input spectrum with an exponential 
cutoff at lO'^ eV. Data points show the atmospheric neu trino background as measured by the AMANDA experim ent (data betwee n 
2000 and 2003) [Munich et all (l2007h : iMiinichl (l2007h . The prediction of atmospheric neutrinos is taken from IVolkoval (Il980l) . 
AMANDA limits are for the same data sample, derived fro m the fact that no signific ant excess above the atmospheric background 
was observed. The dot-dashed line shows the prediction bv iLoeb & WaxmanI ( l2006l) . not taking into account that only 14% of all 
starbursts are calorimeters. 



5.2. Gamma Ray Bursts 



As starburst galaxies show an enhanced rate of supernova ex- 
plosions, an increased rate of long Gamma Ray Bursts (GRBs), 
which are directly linked to SN-Ic events (Mazzali et al., 2003, 
e.g.), is expected. Thus, if long GRBs are the dominant sources 
of UHECRs, the contribution from nearby objects should fol- 



low the distribution of starburst galaxies. Here, we examine the 
number of GRBs to be expected from our catalog. 
In the following calculations, we assume that every SN-Ic ex- 
plosion is accompanied by a particle jet along the former star's 
rotation axis, i.e. by a GRB. The opening angle of the GRB jet 6 
determines, how many SN-Ic can be observed as GRBs, 

"GRB = e ■ «SN-Ic ■ (45) 
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Here, «grb is the GRB rate in a galaxy and e = (1 - cosff) is 
the fraction of SN-Ic producing GRBs. The jet opening angle is 
difficult to determine. Typically, one expects an o pening angle 
of less t han 10° for the prom pt emission, see e. g. iBerger et alj 
(12003b : iRacusin et al.l (l2008l) . Afterg low emission and pre cur- 
sors can have larger opening angles jMorsony et al.L ITOOTI) . As 
we focus on the prompt emission, we will use a typical opening 
angle of ~ 10° as an optimistic estimate, i.e. 



e ^ 0.015. 



(46) 



Further, observational data show that core collapse supernovae 
of type SN-Ib/c contribute with 11% to the total SN rate in star- 
bursts dCappellaro & Turattoll200ll) . Thus, using Equ. (05]), the 
GRB rate in a starburst galaxy is directly correlated to the super- 
nova rate «sn, 

"GRB = e ■ ^ ■ nsN , (47) 

with f ~ 0. 11 as the fraction of heavy SN explosions in all SN 
explosions in a single galaxy. Using Equ. ( |43T l to determine the 
supernova rate in a galaxy in Equ. ( l47T i yields a GRB rate of 



"GRB 



3.8.10--.(^l.(^l.f^l yr 
\ Lo / lo.OlS/ lo.ll/ ^ 



(48) 



per starburst. With an expected lifetime of more than 10 years 
for a neutrino detector like IceCube, luminosities of around 
3 • 10'^ ■ Lq ~ 10'*^ erg/s are required for the detection of a sin- 
gle event. None of the sources in our catalog provides such 
high luminosities: IRAS 17208-0014 is the intrinsically strongest 
source with LpiR = 5.9 ■ lO'*^ erg/s, the second strongest one is 
1RASF17207-0014 with LpiR = 5.5 • lO^^ erg/s. If, however, a 
larger number of starbursts is considered for an analysis, the to- 
tal luminosity increases and with it the probability of observing 
a GRB. Figure [To] shows the total GRB rate for a number of 

A^starbursts galaxicS, 



^Starbursts 

«GRB(^starbiirsts) = ^ nGRB(«th starburst) . 

!=1 



(49) 



In the figure, we sum up the GRB rates achieved in the sin- 
gle starbursts, starting with the most luminous source, adding 
sources in descending luminosity order. The red points show the 
GRB rate summing over all starbursts in the sample. On total, 
0.03 GRBs per year are expected to be observable from our sam- 
ple of starbursts. The blue squares display the total GRB rate, 
summing up sources in the northern hemisphere, which corre- 
sponds to IceCube's Field of View (FoV), in order to estimate the 
neutrino detection probability in Section |5.2.2| Here, 0.02 GRBs 
per year are expected. This number can be enhanced signifi- 
cantly when taking into account those weaker sources which we 
do not include in our catalog in order to ensure completeness 
(see Section[3]l. 

5.2.1 . Observation of cosmic rays and starburst galaxies 

Of the 127 sources in our catalog, 96 are located in the northern 
hemisphere. Thus, Auger South will only have very few star- 
bursts in its FoV, while HiRes, the Telescope Array (TA) and 
later Auger North will be able to observe a large fractio n of 
starbursts. As already discussed by iBiermann et alj (|2008|) . the 
question of the origin of UHECRs can only be resolved by tak- 
ing this strongly asymmetric distribution into account. Even if 
GRBs contribute to the total flux of cosmic rays, nearby sources 
cannot be identified by a southern hemisphere telescope, since 
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Fig. 10. Number of GRBs per year in the starburst catalog, in- 
cluding Nstarimrsts sourccs. Starting with the strongest one. Red 
circles include all starbursts in the sample. The total GRB rate in 
the sample, including all 127 sources, is 0.03 yr"', which means 
that a GRB could be observed every 30 years on average. Results 
shown in blue squares only show northern hemisphere bursts, ly- 
ing in the IceCube FoV. Northern hemisphere GRBs would oc- 
cur every 50 years (;ic«B = 0.02 yr"'). 



the number of sources is too small. Given the large number of 
sources in the northern hemisphere, telescopes like HiRes, TA 
and Auger North can investigate this matter. In particular, for 
GRBs coming from one of the starbursts in our catalog could 
give the opportunity of an enhanced UHECR flux within a short 
time window, assuming that the signal is still focused in time and 
space due to the closeness of the sources. Thus, an analysis with 
UHECR detectors in the northern hemisphere for nearby GRBs 
could be optimized by not only looking for spatial but also for 
temporal clustering. 



5.2.2. Enlianced neutrino flux from GRBs in starbursts 

While the detection of a permanent, diffuse signal from GRBs 
in nearby starbursts will not be possible due to the high atmo- 
spheric background seen by high-energy neutrino telescopes, a 
timing analysis may be able to identify the Gamma Ray Bursts 
in neutrinos. In such an analysis, the location of a nearby star- 
burst can be chosen as a potential neutrino hot-spot. By defining 
a time window of the typical duration of a long GRB (~ 100 s), 
the atmospheric background can be reduced to close to zero. 
The GRB event rate per starburst galaxy was already examined 
above. Here, the general neutrino intensity and in particular the 
possibility of detection with IceCube are discussed. 



Neutrino spectrum from a single GRB in a starburst The neu- 
trino spectrum during the prompt photon emission phase in a 
GRB was for the first time determined by IWaxman & Bahcall 
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([19971 fT999h . It can be expressed as a broken power-law, 



dNy 
dEv 



■ j4y ■ Ey 



for < £v < < (50) 
for By > et . 



The spectrum includes the five parameters Oy, f5y, e'^, el. Ay. 
We discuss the numerical values of th e parameters below. 
For th eir derivation, see e.g. iGuetta et al.l ( |2004|) : iBecker et al.l 
(|2006|) . The derivation is done for an isotropically emitted sig- 
nal. However, the neutrino spectrum itself does not vary with 
the jet's opening angle, since the parameters are determined 
by energy densities. Where the energy is enhanced by a factor 
1/(1 - COS0), the solid angle decreases just as (1 - cos0) and 
the factors cancel. 

- The break energies and 

The neutrino spectrum can be derived when assuming that 
the protons accelerated along the GRB jet interact with the 
ambient synchrotron photon field. Neutrinos are produced 
in TT^-decays as described in Equ. ( |35] |. The first break en- 
ergy represents the energy required to produce the Delta 
resonance in proton-photon scattering. At higher energies, 
Ey > e*, the neutrino spectrum follows the proton spec- 
tral behavior, fiy ~ ap. At lower energies, Ey < e^, scatter- 
ing becomes less effective and the spectrum becomes flatter, 
Oy ~ Op — \. Due to the transformation from the CM sys- 
tem of the proton-photon interaction in the source into the 
observer's frame at Earth, the first break energy depends on 
the shock's boost factor F :- 10^^ ■ F2.5, the observed pho- 
ton break energy of GRBs and the redshift (z ~ for our 
catalog), 



10-** -(1 +zr^ 



F2 
^ 2.5 



e" 



GeV ~ 3 ■ 10^ GeV . (51) 



y,MeV 



We fix the boost factor to F2.5 
break energy to 



-y.MeV ■' 



e°/MeV 



1 and the observed photon 
0.25. For the sources in 



our catalog, we have z ~Q- 

The second break in the neutrino spectrum, due to pion- 
synchrotron losses, is determined by unknown parameters 
like variability time scale f,, x 0.01 s ■ty-2, GRB luminosity 



10"' erg/s 



tition fractions, e^. 



and electron and magnetic field equipar- 

6B *0.1. 



3 ■ 10^ 
1 +z 



pl/2 -l/2p4 . , 



L^i GeV 



3 ■ 10' GeV. 



the parameters to f,,,-2 = 1, 



(52) 

L^: = 1, 



Here, we fix 

= 6^=0.1. 
The spectral indices Uy and fiy. 

In the energy range < Ey < e^, the spectral index = p 
is represented by the primary hadron spectral index. We as- 
sume a simple E^^ proton spectrum in this case; While it was 
argued previously, that highly relativ istic flows have a hmit 
of no t becoming flatter than ~ E^^-^ ( iBednarz & Ostrowskil 
Il998h . recent studies of oblique shocks have shown that par- 
ticle spectra as flat as Zs"'-^ can be produced with s hock boost 
factors of F > 100 as they occur in GRBs, see iMeli et aT] 
(I2008h . In addition, including large angle scattering of par- 
ticles in the acceleration process also leads to flat particle 
spectra, see IStecker eTa l. (2007). This implies, the exact 
spectral index depends on the shock properties. As we are 
not dealing with observed GRBs here, but with potentially 



to-be-observed GRBs, we cannot determine the spectral be- 
havior by the observation of the synchrotron spectrum as it 
was done injOuetta et al. (2004); Becker et al. (2006). Thus, 
we use Uy - I and jSy = 2 as a first order approximation, just 
as it is done in s everal approaches of unknown sin gle-source 
spectra, see e.g. (IWaxman & Bahcal]Lll997lll999h . 
The normalization Ay 

The normalization, which corresponds to the intensity of the 
burst, is given by 



1 1 



Av= -- 



piso 



fn 



^fe4n-dj In(lO) 



(53) 



Here, 



of fgo 



we 



use a fixed isotropic energy release of 
y ■ tgo - 10^^ erg, where we assume a burst duration 
10 s. Further, we assume that the energy contained 
in GRB electrons is 1/lOth of the energy contained in GRB 
protons, fe - OA. The factor ^ 0.2 describes the fraction 
of energy transferred to the charged pion. The factor 1/8 is 
applied since 1 /2 of all proton-photon interactions go into 
neutrino production, and 1/4 goes into a single neutrino fla- 
vor Here, we fix all parameters except for the distance di of 
the GRB from Earth, which is determined by the distance of 
the starburst from Earth. 

To estimate the neutrino flux from a standard GRB from one 
of the starburst galaxies in our sample, we simply calcu- 
late a normalization, dependent on the distance of the star- 
burst. Every other parameter is kept constant. Hence, the re- 
sults can only serve as a rough estimate of what would hap- 
pen on average. Both the spectral indices, the break ener- 
gies and the normali zation fluctuate with each burst as de- 
scribed theoretically in lGuetta etal] (|2004 : [Becker et al.1 (l2006h 
and worked out experimentally i n AM ANDA for the case 
of GRB030329 ("Stamatiko s et al.l [2005b and in IceCube for 
GRB080319B (Kappes et al.L 120081) . One of the most impor- 
tant parameters is the total GRB energy. Here, we assume that 
it is E"° = 10^^ erg. The isotropic equivalent energy for the 
'naked eye' GRB0803 19B is two orders of magnitude higher, i.e. 
1 .3 ■ 10'''' erg and several weak bursts will have smaller isotropic 
equivalent energies. While these eff'ects cannot be taken into ac- 
count due to lack of knowledge, they should be kept in mind. 
In case of a positive detection in a high-energy neutrino detec- 
tor, the parameters can of course be determined explicitly. In the 
case of a negative result, it would be possible to set a limit in the 
following sense: no GRB with isotropic equivalent luminosities 
occurred during the time of observation in the given starburst(s). 



Expected events in IceCube The probability of a GRB to be 
observed from one of the starburst galaxies in the sample from 
the northern hemisphere is shown in Fig. [T0| (blue squares). 
Analyzing all 96 sources in the northern hemisphere, a rate of 



,. total/noith 

''grb 



(A^sta 



starbursts 



96) = 0.02 yr" 



(54) 



is expected. This number can be enhanced significantly when we 
consider all starbursts originally selected (309 sources). We only 
consider the brightest ones for completeness reasons. Thus, we 
expect that at least one GRB from a starburst in the supergalac- 
tic plane should happen in the lifetime of IceCube, if not more. 
Prospects for Km3NeT are not as optimal as for IceCube, since 
Km3NeT will be located in the northern hemisphere, looking 
at the southern sky, and the instrument will therefore only see 
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Fig. 11. Neutrino intensity Ay for a single GRB from a starburst 
at a luminosity distance di. A generic GRB isotropic energy 
of - 10^^ erg is assumed. Therefore, the neutrino intensity 
only depends on the luminosity distance as A,, oc d^^. 
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Fig. 12. Number of GRBs in a starburst versus neutrino flux in- 
tensity for each starburst. The sharp cutoff" at low neutrino inten- 
sities comes from the distance cut performed at z = 0.03. This 
leads to a cutoff in the neutrino intensity, since the latter is pro- 
portional to the luminosity distance squared. 



Fig. 13. Histogram of number of events in IceCube from the 
96 starburst galaxies in the northern hemispher e. The IceCube 
effective area we use jMontaruli et all l2007h does not take 
into account sensitivity changes with declination. All bursts are 
stronger than 1 event in IceCube, since they come from star- 
bursts closer than z - 0.03. For more distant bursts, the number 
of events will be less, since the neutrino intensity decreases with 
the distance squared. 



Here, we neglect the weak dependence of the effective area on 
the declination of the burst. For th e threshold energy, we use 
£,/, = 100 GeV ( Ahre ns et alll2004l) . This is the general detec- 
tion threshold of IceCube. As events can be selected by direc- 
tion and a small time window, the atmospheric background can 
be reduced to almost zero. Therefore, events in the entire energy 
range are available in such an analysis. If we now assume a stan- 
dard burst with an isotropic energy release of Ey° - 10^^ erg, 
we can estimate how strong a GRB would be on average from 
the nearby starbursts in our catalog. Figure [13] shows the his- 
togram of the number of events per starburst for an average 
GRB. Depending on the distance of the starburst, the event rate 
in IceCube ranges from more than 1 event up to more than 1000 
events per burst in a small time window of around 10 - 100 sec- 
onds. These numbers lie between 2 and 6 orders of magnitude 
above those GRBs typically observed by satellite experiments 
like BATSE, Swift and Fermi. The reason why such a strong 
burst has not been observed in gamma-rays is simply the rela- 
tively low rate of occurrenc^B If such an event happens, neutrino 
detectors, and also cosmic ray detectors have a higher chance of 
detection due to their extremely large FoV (~ 2n). 



a small fraction of the class of starburst galaxies, which domi- 
nantly shows sources in the northern hemisphere. 
IceCube's effective area A^/f is presented by iMontaruli et alj 
(l2007h and can be used to determine the total number of events 
per GRB in IceCube, Nevems, by folding it with the GRB spec- 
trum dNv/dEy, 

Ne.en,s^ f A,ff(Ey)-^{Ey)dEy. (55) 

Je,„ fl^v 



5.2.3. High-energy photons and nearby GRBs 

Following the channel of tt" production rather than the neutrino- 
producing charged pions, we see that a high-energy photon sig- 
nal is always expected along with high-energy neutrino emis- 
sion. Such a correlation could be observed for sources optically 



^ The exact number cannot be given here, since we only consider 
those sources with relatively high fluxes in the radio and at FIR wave- 
lengths. The total rate of 0.02 GRBs/year could be enhanced signifi- 
cantly with a lower flux threshold. 
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thin to TeV photons. In such a case, GRBs in nearby starburst 
galaxies would produce an equally strong signal in neutrinos as 
in high-energy photons. This means a good discovery potential 
for large FoV detectors at GeV-TeV energies. Fermi is one of the 
best candidates to observe this high-energy radiation at Ge V en- 
ergies, Milagro and its successor HAWC will be able to see such 
events at TeV energies. Air Cherenkov Telescopes like H.E.S.S., 
MAGIC and VERITAS have a FoV that is too small and are 
therefore not well-suited for a detection of those rare GRBs. If, 
however, such a GRB is detected by Swift, especially MAGIC 
will be able to follow up the GRB quickly, as the telescope is de - 
signed for quick GRB follow-ups, see e.g 



5.2.4. Summary 



as the telescope is de - 
iGalanteet al.1(l2008l) . 



To sum up, if a GRB will happen in one of the 127 starburst 
galaxies, the detection probability is extremely high in wide- 
FoV detectors for high-energy neutrinos and photons as well as 
in cosmic rays. Since most sources are located in the northern 
hemisphere, detectors like IceCube, HAWC and Auger are opti- 
mal for for such a study. This would enable a detailed study of 
the hadronic emission processes of GRBs, since the main three 
messenger particles are covered. Explicit search methods for an 
enhancement of cosmic rays, high-energy photons and neutrinos 
from those starbursts can clarify weather or not such a burst hap- 
pens or not. In the case of the closest sources, the GRB should 
even be visible in a typical skymap of neutrino, photon or cosmic 
ray telescopes. 



6. Implications and possible experimental tests 

In this paper, we present a model to explain the correlation be- 
tween far-infrared and radio emission in starburst galaxies and 
we discuss the particle spectra from cosmic ray interactions and 
possibilities of their observation at Earth. In this context, we also 
present a catalog of 127 nearby (z < 0.03) starburst galaxies with 
both FIR (560^ > 4 Jy) and radio (S i.4ghz > 20 mJy) measure- 
ments. Those starbursts lie in the supergalactic plane (z < 0.03) 
and serve as a test for basic properties of starbursts. It can further 
be used to calculate the cosmic ray emission from starbursts in 
the local Universe. In this section, we summarize implications; 

I. FIR-Radio correlation 

Within this model, we show that the radio emission is basi- 
cally independent of the magnetic field strength. This leads 
to a calorimetric correlation between radio and FIR emis- 
sion. 

Starburst galaxies as well as quiescent galaxies equally obey 
the correlation between the radio and the far-infrared emis- 
sion. Therefore the theory implies by necessity that the mag- 
netic field in a starburst region rises with the overall energy 
density. Since galactic dynamos are not understood, this just 
raises the ante in terms of the requirements for any the- 
ory that tries to account for the magnetic fields and their 
strengths in galaxies, see ( iBeck&HoerneslI 19961 e.g.). 
Interacting galaxies, subject to substantial compression, 
show magnetic fields higher than corresponds to the energy 
density derived from star formation. Therefore, it may be 
expected that they have more radio emission than indicated 
by the general radio/FIR cor relation. At least one example 
shows this kind of behavior iH ummgl & Beck, 1995). 
There is likely also a lowest level of star formation activity, 
where the assumption that the hot medium is fully connected. 



fails. There one may expect also substantial departures from 
the correlation. 

High-energy photons and neutrinos from SNRs in starbursts 
The diffuse high-energy neutrino and photon flux from 
SNRs in starbursts h a s been discussed pre viously by 
iLoeb & Waxmanl(l2006l) : lThompson et al.l(l2007l) . In our cal- 
culation, we assume that a fraction of the total energy of a 
SNR is transferred to cosmic ray acceleration. It is usually 
assumed that supernova remnants in the Galaxy accelerate 
particles up to the knee, i.e. Zsmax ~ 10'^ eV. Under the as- 
sumption that this is the maximum energy for SNRs in star- 
burst galaxies, we show that a detection of the diffuse neu- 
trino flux from starbursts with IceCube or Km3NeT is un- 
likely. The main reason is the low energy cutoff, and not the 
strength of the signal. On the other hand, if it were possible 
to accelerate protons to e ven higher energies in SN Rs in star- 
bursts, as is suggested bv iLoeb & WaxmanI (|2006|) . IceCube 
will be sensitive to the neutrino flux. High-energy photons 
make up around 10% of the total high-energy photon back- 
ground observed by EGRET. The high-energy photon de- 
tector Fermi was launched on June 1 1 , 2008 and will soon 
be able to give more details on the background and which 
sources and source classes contribute. Cosmic rays them- 
selves cannot be observed from starburst galaxies, since the 
background of Galactic cosmic rays is too high. 
High-energy neutrinos and cosmic rays from GRBs in nearby 
starbursts 

Due to an enhanced rate of star formation and supernova ex- 
plosions in starbursts, these Galaxies provide a higher prob- 
ability for Gamma Ray Bursts compared to regular galaxies. 
So far. Gamma Ray Bursts could not be identified in high- 
energy neutrinos. Despite their tremendous release of neutri- 
nos, they are usually too far away to yield a significant sig- 
nal individually (z ~ 1-2). The stacking of known GRB 
locations in ne utrino data was done with the AMANDA 
experiment, see lAchterberg et al.l (l2007i l2008l) . with limits 
above the typical, expected diffuse flux. IceCube will be 
able to explore the diffuse neutrino flux using the stacking 
method. Two single GRBs were also anal ysed, the first one 
being GRB030329, using AMANDA data dstamatikos et al.L 
I2005i) . the second one bein g GRB080319B, us ing IceCube 
in a 9-string configuration dKappes et al.L l2008h . In this pa- 
per, we show using our catalog of starbursts in the super- 
galactic plane, that GRBs in closeby starbursts (z < 0.03) 
typically produce event rates in IceCube above 1 event per 
bursts. For the most distant sources at z = 0.03, about 1 
event is expected, for the closest sources, more than 1000 
events are predicted within a short time window of a few 
tens of seconds. In our catalog of 127 starbursts, about 1 
GRB is expected each 30th year, most of those happen- 
ing in the northern hemisphere. The total rate of GRBs in 
the supergalactic plane is, however, expected to be higher, 
since we only include the brightest sources in this catalog 
(560// > 4 Jy and 5i.4ghz > 20 mJy). The stacking of 
starburst galaxies for a neutrino analysis would help iden- 
tifying the sources at z ~ 0.03, while the closest sources 
will immediately be visible in a neutrino skymap. Those 
GRBs should be accompanied by high-energy photon emis- 
sion as well. Fermi will be able to identify these at GeV ener- 
gies ( Connaughtoni 2007h. At even hig her energies, Milagro 
JVasileiou & (MILAGRO Coll.l |2007[) and th e next gener- 
ation experiment HAWC dSinnis et al.L 120051: ISmith et all 
'2006*) will be able to investigate TeV emission. For such 
nearby and strong cosmic ray emitters, it might even be pos- 
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sible to identify them in charged cosmic rays, given that dif- 
fusion in time and space is not too large. 
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Appendix A: Catalog tables 

We present the catalog of nearby starburst galaxies discussed 
throughout the paper The 127 sources presented here were se- 
lected from a larger sample of 309 sources, all previously iden- 
tified as starburst galaxies. We require the FIR flux at 60 fi to 
be larger than S(,q^i > 4 Jy and the radio flux at 1.4 GHz to 
be larger than S iaghz > 20 mJy. In addition, we only include 
sources at z < 0.03. Table lA. ll summarizes the basic properties of 
the catalog: name, right ascension (RA [deg]), declination (DEC 
[deg]), redshift (z) and distance (Dl [Gpc]), together with the 
FIR flux measurements, S a, A giving the wavelength. Table |A.2l 
presents radio flux measurements, Sy, with v as the frequency. 
Table lA. 3l summarizes X-ray flux measurements. References are 
always given in the last column. 



Table A.l. Coordinates, distances and IRAS measurements of the sample. All fluxes in [Jy]. Coordinates epoch J2000.0 and distances are obtain ed from NED. In NED, distances are corrected to 
the cosmological microwave baclcground using cosmological parameters according to ACDM cosmology, ho = 0.73, Q.,,, = 0.27, Da = 0.73, using lWrightl l l2006h . 

I RAS Reference s: 

l: ISandersetaD l l2003h . 2: lMoshir et aTHl990bh . 3: lLisenfeld etaP l l2007h . 4: ISurace et all ( l2004h . 5: lKnappl l fT994h . 6: lBeichman et alj ( ll988h . 7: ISoifer et an ( ll989h 



Name 


RA [deg] 


DEC [deg] 


z 


mCpc] 


Sl2/jm 


S25/jm 


S60/jm 


SlOO/jm 


References 


MRK545 


2.47254 


25.9238 


01 523 


05962 

\j .\j ^ y yj ^ 


0.523 


1.082 


9.196 


15.34 


2 


NGC34 


2.77729 


-12.1073 


019617 

\j .\j 1. yyj 1. 1 


0.0771 


0.35 


2.39 


17.05 


16.86 


1 


MCG-02-0 1-051 


4.71202 


-10 3768 


027103 

yj m yj ^ 1 1. yj ^ 


0.109 


0.24 


1.19 


7.35 


10.22 


4 


NGC174 


9.24558 


-29.4778 


01 1905 


0.0451 


0.41 


1.27 


11.36 


19.77 




NGr232 


10.6909 


-23.5614 


022172 


0886 

\j .\jyjyjyj 


0.36 


1.28 


10.05 


17.14 




NGC253 


11.888 


-25.2882 


0008 

yj ^yjyjyjyj 


0.0031 


41.04 


154.67 


967.81 


1288.15 




IC1623 


16.9466 


-17.507 


02007 


07857 


1.03 


3.65 


22.93 


31.55 




NGC520 


21.1461 


3.79242 


00761 

yj *yjyj i yj ±. 


03022 

vJ ■ v ^ yj ^ ^ 


0.9 


3.22 


31.52 


47.37 




NGC632 


24.323 


5.87764 


010567 

yj ■ \j ^ \j -y yj i 


0.0396 


0.37 


0.88 


4.89 


7.32 




NGr660 


25.7598 


13.6457 


00283 

\J .\J\J ^yj ^ 


01233 


3.05 


7.3 


65.52 


114.74 




NGr828 


32.5399 


39.1904 


0.01793 


07073 


0.72 


1.07 


11.46 


25.33 




NGC891 


35 6392 


42.3491 


00176 

yj *yjyj ±. i yj 


00857 

yj »yjyjyj^ i 


5.27 


7 


66.46 


172.23 




NGr958 


37 6785 


-2.939 


0.01914 


0.0765 


0.62 


0.94 


5.85 


15.08 




NGn055 


40 4385 


443 1 67 


00332 


0.01131 


2.24 


2.84 


23.37 


65.26 




Maffei2 


40.4795 


59.6041 


-5.7e- 05 


00332 

yj *yjyj ^ ^ ^ 


3.624 


9.238 


135 


225 




NGri068('M77~) 


40 6696 


-0 0132806 


00379 


0.0137 


39.84 


87.57 


196.37 


257.37 




UGr2238 


41.5729 


13.0957 


021883 

\J .\J^ -I yjyj ^ 


0883 


0.36 


0.65 


8.17 


15.67 




NGn097 


41.5794 


-30.2749 


00424 


0.0152 


2.96 


7.3 


53.35 


104.79 




NGCl 1 34 

± 1 v_> J. J. ^ r 


43.4222 


13.0141 


012142 


0.0474 


0.55 


0.92 


9.09 


19.43 




NGCl 365 

J- 'I V_> J. ^ vj»y 


53 401 5 


-36.1404 


00546 


0.01793 


5.12 


14.28 


64.31 


165.67 




IC342 


56.7021 


68 0961 


0.0001 


0.0046 


14.92 


34.48 


180.8 


391.66 




UGr02982 


63 0935 


5.54739 


017696 

\j .\j 1. 1 yj y yj 


0.0724 


0.57 


0.83 


8.39 


16.82 




NGCl 530 

1 > V.J J. »y 


65 8629 


75.2956 


00821 


03622 


0.72 


1.23 


9.88 


25.88 




NGCl 569 


67.7044 


64.8479 


-0 00035 

yj ^yj yj yj ^ ^ 


0.0046 


1.24 


9.03 


54.36 


55.29 




MRK617 


68.4994 


-8 57888 


0.01594 


06261 


0.441 


7.286 


32.31 


32.69 




NGCl 672 

± 1 V_ J V.^ J. \J / ^ 


71.4271 


-59.2473 


0.00444 


01682 


2.47 


5.25 


41.21 


77.92 




MRK1088 


73 6598 


3.26797 


01528 

yj *yj 1- ^ ^ yj 


0605 1 

yj *yjyjyj^ ±. 


0.2659 


0.835 


6.605 


10.77 




NGCl 808 

1 1 V.J ^ — . J. u v/u 


76.9264 


-37.5131 


00332 


01261 


5.4 


17 


105.55 


141.76 




NGCl 797 


76.937 


-8 01908 

yj .\j 1. y\jyj 


014814 

\j m \j 1. r yj 1. r 


0.0616 


0.33 


1.35 


9.56 


12.76 




MRK1194 


77.9423 


5 20061 

^ t^yj y^yj ±. 


0.01491 


05948 

yj ^yj -y y ^yj 


0.283 


0.7071 


6.688 


11.5 




NGC2 1 46 


94.6571 


78.357 


00298 

\j .yjyj ^ y yj 


0.012 


6.83 


18.81 


146.69 


194.05 




NGC2276 


111.81 


85 7546 

/ ^ ryj 


00804 

yj *yjyjyj yj r 


0.0328 


1.07 


1.63 


14.29 


28.97 




NGC2403 


114.214 


65.6026 


0.00044 


0.00247 


2.82 


3.57 


41.47 


99.13 




NGC2415 


114.236 


35.242 


0.01262 


0.05341 


0.61 


1.19 


8.75 


13.58 




NGC2782 


138.521 


40.1137 


0.00848 


0.03951 


0.64 


1.51 


9.17 


13.76 




NGC2785 


138.814 


40.9175 


0.008746 


0.0392 


0.49 


1.09 


8.4 


15.79 




NGC2798 


139.346 


41.9997 


0.00576 


0.02784 


0.76 


3.21 


20.6 


29.69 




NGC2903 


143.042 


21.5008 


0.00186 


0.00826 


5.29 


8.64 


60.54 


130.43 




MRK708 


145.548 


4.67314 


0.00682 


0.03116 


0.46 


0.8 


5.36 


8.24 




NGC3034(M82) 


148.968 


69.6797 


0.00068 


0.00363 


79.43 


332.63 


1480.42 


1373.69 




NGC3079 
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Table A.2. Radio measurements of the sample. All fluxes in [mjy]. 
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_ 


1,2,8,14 


NGC34 


67.5 


— 


_ 


— 


— 


_ 


_ 


4 


MCG-02-0 1-051 


43.2 


— 


_ 




_ 


_ 


_ 


4 


NGC174 


45.7 


_ 


_ 


_ 


_ 


_ 


_ 


4 


NGC232 


60.6 


— 


— 


— 


56 


— 


— 


4 


NGC253 


6000 


_ 


— 


3520 


2433 


— 


2580 


5,16,17 


IC1623 


249.2 


_ 


— 


— 


96 


_ 


_ 


4,5 


NGC520 


176 


110 


_ 


— 


87 


_ 


_ 


1,2,8 


NGC632 


23 


15 


— 


— 


— 


_ 


_ 


2,8 


NGC660 


373 


255 


_ 


_ 


187 




_ 


1,4,8 


NGC828 


108 


_ 


_ 


— 


47 


_ 


_ 


6 


NGC891 


701 


— 


— 


— 


342 


— 


— 


2 


NGC958 


71.9 


_ 


_ 


_ 


_ 


_ 


_ 


4 


NGC1055 


200.9 


129 


_ 


150 


63 


_ 


_ 


1,4,8,16 


Mafrei2 


1015 


_ 


_ 


— 


375 


_ 


— 


1,19 


NGC1068(M77) 


4850 


_ 


_ 


305 


2039 


1890 


1342.4 


2,9,13,14,17 


UGC2238 


72.2 


— 


— 


— 


— 


— 


— 


1,2 


NGC1097 


415 


— 


— 


250 


126 


— 


150 


6,7,16 


NGC1134 


89.1 


57 


_ 


_ 


32 


_ 


_ 


1,2,8 


NGC1365 


530 


_ 


_ 


350 


230 


180 


— 


4,7,16 


IC342 


2250 


_ 


_ 


— 


277 


_ 


— 


1,6 


UGC02982 


91.3 


_ 


_ 


_ 


_ 


_ 


_ 


2 


NGC1530 


80.7 


— 


— 


— 


27 


— 


— 


1,6 


NGC1569 


396 


— 


— 


— 


198 


— 


155 


1,10,19 


MRK617 


138 


_ 


_ 


_ 


63 


_ 


_ 


4,13 


NGC1672 


450 


_ 


_ 


210 


114 


_ 


100 


16,18 


MRK1088 


45.7 


31 


_ 


_ 


_ 


_ 


_ 


2,8 


NGC1808 


497 


_ 


_ 


350 


229 


_ 


220 


6,16,18 


NGC1797 


29.1 


- 


- 


- 


- 


- 


- 


4 


MRK1194 


42.2 


27 












2,8 


NGC2146 


1087 














6 


NGC2276 


283 














6 


NGC2403 


387 








169 






19 


NGC2415 


66.4 


53 






41 




30 


1,2,8,10 


NGC2782 


107.5 








47 




60 


1,10,20 


NGC2785 


67.6 














2 


NGC2798 


82.0 








37 




53 


1,2,10 


NGC2903 


444 


200 






118 






1,2,8,13 


MRK708 


32.6 


21 












2,8 



CO 



n 

o 



73 



Table A.2. continued. 



Name 


Si dnCHy 

J. . VvXlA 




Si liQnHy 


Si 70 HHy 


Sd CH? 


S^ no ClWr 


S^ ni CHy 


References 




/ZoD.O 












■^Ql 7 
1 z 


1 1 7 90 

1,1/, Zu 




870 7 








^91 
Jz i 






1 91 

1 , Z 1 












AA 




8 1 

0. 1 


1 9 9^ 

1 , Z, ZJ) 














740 
Z^W 


7S0 

ZJW 


J, O, 1^, lO 


MRK33 


24.6 


















417 








1 59 






1 1 9 


i> vJ^J JO / 


1 1 8 


71 












1 9 8 1 1 

1 , Z, o, lU, ID 
















jy 


1 1 
1, IW 




11 A. 

Z / T- 








1 1 7 






1 9 8 
1 , z, o 




Z^J 








76 
/ D 






1 1 Q 

1, IV 






900 






141 
1^1 






1 9 8 
1 , z, o 




/L7n ? 

T- / W.Z 








776 

z / VJ 


700 
zuw 


774 

ZZH- 


1 8 10 1 /L 9 1 

1,0, IW, lT-,Zl 




1 97 

iZ / 














z 
















"^67 

JDZ 


1 9^ 

lU, Zj 


iviivrv i o o 


^0 7 












7S 


9 1 

Z, 1 w 




ijy 














1 9 
1 , z 




70 8 


SO 






S7 

JZ 






1 9 8 
1 , z, o 




1/L7 

It- / 










S6 




1 zL \f, 1 Q 

IH, lU, 1" 












48 






1 9 

i, z 


MnPdl 09 
i> vJv^T- i UZ 


Z. I J 








70 
/ U 




1 OS 
1 WJ 


1 1 

1, lU 


MR K 1466 


90 


1 5 












2 8 

z, o 


iVlKJV 1 Jy 


"^1 Q 


1 fs 






1 

1 J 






9 8 1 S 

Z, o, 1 J 




iZZ 














1 1 Q 

1, IV 


NHPzL? 1 A 


^8 ^ 
J>o . J) 














9 
z 


IN vJv^'t-Z / J 


78 S 


DJ 






"^7 
J / 






9 8 1 S 

Z, o, 1 J 




/I /I /I 








1 07 
lUZ 


1 70 
IZW 




9 8 IzL 1 S 

Z, O, It-, i-J 




777 
zz / 


1 ^8 






7S 

/ J 






1 8 1 q 

1,0, 1" 




"^8 S 














9^ 




1 87 Q 


1 7Q 






77 






9 8 1 S 

z, o , u 


NnP4S^6 


704 Q 








1 1 4 

1 IH- 


1 1 
1 1 w 




9 8 1^ 

Z, O, 1 J 




1 1 77 
i i zz 








4'^8 






1 8 1 Q 
1, o, IV 


NnPAf^66 


A'KA 








161 

lO 1 






9 1 ^ 

Z, 1J> 




1 1 ^ 


72 






46 






12 8 


^C^CASKlf\(\Af\A\ 
i> vJv^T-OZDViViO'T- ^ 




D / 






S4 






1 9 8 
1 , z, o 




DOUW 






soon 






7840 


^f^ 18 

ID, lo 




1 Q4 








67 

VJZ 






1 1 q 

1,1V 






77 
zz 












9 8 
Z, 6 












1 74 






1 9 

i, z 


ARP193 


104 








53 






1 2 

J. , 


NGC5104 


39.9 


38 












2,8 


NGC5135 


194 








107 






6,7 


NGC5194(M51) 


1310 








436 




360 


1,10,19 


NGC5218 


30.4 














2 


NGC5236(M83) 


2445 








648 






6,7,16 



Table A.2. continued. 



Name 


Sj 40 GHz 


S2.38 GHz 


S2.69 GHz 


S2.70 GHz 


S4 85 GHz 


S5 00 GHz 


S5 01 GHz 


References 












Al 






i, ly 






AS 












7 8 
z, 




O J .0 








QO 

y\j 


7S 




7 1/1 


1 Tr;p87^Q 

Uvj^o / jy 


yj .u 








J / 






1 9 
1 , Z 


MR 1^1 ^fiS 


.yj 


1 A 












9 8 
Z, 












9Q 

Ly 




/LO 


i , Z, i u 


















Of, 
ZD 


i>vJ^JD / O 


1 1 n 








f^8 
Do 






z 




1 1 f, 

i iU 








^8 

JO 






9 
z 




997 

zzz 








V J 


1 J 




i , Z, 




991 

ZZ i 


1 "^8 






D / 






1 8 1 Q 




UU.H- 










1 7 




9 1 

Z, i w 






81 
i 






ASK 
H-O 






19 8 10 

i, Z, 0, iU 


ARP9?n 
rVlvr ZZW 


JZ'4- 




9^0 
ZDU 




908 
ZWO 






1 9 ^ 

1 , Z, J 


i> vJv^J"UZ 


89 ^ 

OZ. J 


JO 






JO 






1 7 8 
1 , z, 




f>'^ Q 


oy 












9 8 
Z, 


Mnpfi 1 8 1 

1 > vJ^D i O i 


yj .u 








JO 






1 9 8 
i , Z, 


i> vJv^UZ i / 


7Q Q 
ly .y 














9 
z 


i>vJ\^DZM-u 










1 7Q 

i ly 




1 70 


1 ^ 1 1 Q 


INvJ^DZou 


i J / 








J J 






1 9 
1 , z 


ll\./\o iOZ7J J)T-U 


99^ 

ZZ J) 








1 AA 






7 


Nr;pfi7ni 


Q9 9 
yz.z 








99 

ZZ 








Nr;pfi7f^zL 


1 1 s n 








'XA 
J^ 






19 10 
i, Z, iu 












1 

J J i 






1 1 q 


MripTi an 

INvJ^ / i J\J 


io J 














/; 



TPS17Q 


iUJ 








7Q 

/ y 






fi 7 


wr;p7^^i 

1 > vj^ / J) J) i 


J / J 


1 87 

10 / 






80 

ou 






1 9 8 
1 , Z, 


Nr;P7zLfiQ 


9SS 


1 






VJ 






8 1^10 
0, IJ), Ly 


Nr;P7zl7Q 


107 


S8 

Jo 






zll 






1 9 8 
1 , Z, 


















fi. 


\rr;r'7S/Li 


lf^9 
iOZ 


1 01 






JO 






1 9 8 

1 , Z, 


i^JZVO 


9ii 9 

Z'+.Z 














9^ 

Z J 




97f^ 
z / u 








1 

i J" 






1 8 
0, 1 


i> vj^ / Jy 1 


S9 1 

JZ. i 


JV 












9 8 
z, 


Nr;p7SQ9 

i> vj^ / Jyz. 


7S 

/ J 














97 

z / 


MRK^^ 1 Q 


J i .D 


91 

Z i 












9 8 
Z, 


ly V Jv , / D / 




"^0 
JW 












9 8 
Z, 


1 > V.J v.. / D / o 


^y .J 


JO 












9 8 

Z, 


MRK534 


55.8 


33 






45 






2,8, 13 


NGC7679 


55.8 


33 






45 






2,8,13 


NGC7714 


65.8 








39 






1,2 


NGC7771 


229 


97 






57 






1,8,19 


NGC7793 


103 














6 


MRK332 


36.5 


27 












2,8 



Table A.3. X-Ray measurements of the sample. All fluxes in [nJy] . 
References: 

1: Rigopoulou et all ( Il996h. 2: iFabbiano et alJ ( Il992h . 3: IWhite et alUlOOOh . 4: iTaier et alj ( l2005h . 5: lBrinkmann et alj ^1994^ , 6: lOtt et alj ( |2005|) , 7: iTiillmann etai] ( l2006h , 8: lTeng etal.U2005h . 
9: fShu etaLl ( l2007h . lo TGuainazzi et al.l ( l2005h 



Name 


EO IPC(ROSAT) 


EINSTEIN 


Chandra 


ROSAT 


ROSAT 


Chandra 


XMM 


References 




0.1-4 keV 


0.2-4.0 keV 


0.1-2.4 keV 


0.2-2.0 keV 


0.1-2.4 keV 


0.3-8 keV 


0.3-2 keV 




NGC34 














23 


9 


NGC520 








33.4 








3 


NGC660 








45.1 








3 


NGC891 














69.4 


7 


NGC1068(M77) 




3940 






10900 






2,5 


NGC1097 




578 












2 


NGC1365 




326 












2 


IC342 




939 












2 


NGC1569 




368 












2 


MRK617 




124 












2 


MRK1088 








42.3 








8 


NGC1808 










446 






5 


NGC2403 




364 












2 


NGC2415 








63.8 








3 


NGC2903 




273 




46.4 








2,3 


NGC3034(M82) 




4490 












2 


NGC3079 




113 












2 


NGC3256 










1060 






5 


NGC3310 




208 












2 


NGC3367 


102 














1 


NGC3448 




72.6 












2 


NGC3504 




76.8 












2 


NGC3627(M66) 










603 






5 


NGC3690 




91.5 












2 


NGC4102 














45.1 


9 


NGC4214 












75.6 




6 


NGC4273 








67.1 








3 


NGC4303(M62) 




180 






233 






2,5 


NGC4536 




112 












2 


NGC4631 




263 












2 


NGC4826(M64) 




150 




53.3 








3 


NGC4945 












794 








5 


NGC5135 




63.9 




176 








3 


NGC5236(M83) 




972 






954 






2,5 


NGC5256 














9.93 


10 


NGC5253 




41.4 




38.5 




35.8 




2,3,6 


NGC5775 














15.8 


7 


ARP220 


45.69 














1 


NGC6240 










914 






5 


NGC6946 




611 












2 



CO 



n 
o 



73 



Table A.3. continued. 



Name 


EO IPC(ROSAT) 


EINSTEIN 


Chandra 


ROSAT 


ROSAT 


Chandra 


XMM 


References 




0.1-4 keV 


0.2-4.0 keV 


0.1-2.4 keV 


0.2-2.0 keV 


0.1-2.4 keV 


0.3-8 keV 


0.3-2 keV 




i> V J V. / JJ L 




iu.o 












Z, J 






1 7000 














NGC7552 




167 












2 


MRK534 




204 












3 


NGC7679 




204 












2 


NGC7714 




53.8 












2 


NGC7771 




138 












2 


NGC7793 




130 












2 



S' 
CO 



n 

o 



73 
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